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■ ABSTRACT 

o ■ 

Both the nonrelativistic scattering and the spectrum in the presence of the Aharonov- 
Bohm potential are analyzed and the single-particle density of states for different self- 
^ \ adjoint extensions is calculated. The single-particle density of states is shown to be a 

\ symmetric and periodic function of the flux which depends only on the distance from 

the nearest integer. The Krein-Friedel formula for this long-ranged potential is shown to 
be valid when regularized with the C, function. The limit when the radius R of the flux 
tube shrinks to zero is discussed. For i? 7^ and in the case of an anomalous magnetic 
moment Qm > '2, (note, e. g., that Qm = 2.002 32 for the electron) the coupling for spin- 
down electrons is enhanced and bound states occur in the spectrum. Their number does 
depend on a regularization and generically does not match with the number of zero modes 
in a given field that occur when = '2. Provided the coupling with the interior of the flux 
tube is not renormalized to a critical one neither bound states nor zero modes survive the 
limit i? — > 0. The Aharonov-Casher theorem on the number of zero modes is corrected 
for the singular field configuration. Whenever a bound state does survive the R —>■ 
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limit it is always accompanied by a resonance. The presence of a bound state manifests 
itself in the asymmetric differential scattering cross section that can give rise to the Hall 
effect. The Hall resistivity is calculated in the dilute vortex limit. The magnetic moment 
coupling and not the spin is shown to be the primary source for the phase-shift flip that 
may occur even in its absence. The total energy of the system consisting of particles 
and field is discussed. An application to persistent currents in the plane for both spinless 
and spin one-half fermions is given. Persistent currents are also predicted to exist in the 
field of a cosmic string. The 2nd virial coefficient of anyons with a short range (5-function 
interaction is calculated. The coefficient is shown to be remarkably stable when such an 
interaction is switched on. Several suggestions for new experiments are given. 

PACS numbers : 03.65.Bz, 03-80.-hr, 05.30.-d, 73.50.-h 



(Phys. Rev. A 53, 669-694 (1996)) 



1 Introduction 



In this paper we shall report on several physical phenomena and calculate several quanti- 
ties in the presence of the Aharonov-Bohm (AB) potential which, in the radial gauge, 
is given by 

$ OL 

^ = 0, A^ = — = — $0. (1) 
ivrr zvrr 

Usually, $ = a $0 is the total flux through the flux tube and $o is the flux quantum, 
$0 = /ic/|e|. The AB potential will be considered here in a more general sense since 
formally the same potential (of nonmagnetic origin) is generated around a cosmic string. 
The parameter $ is then l/Quiggsi Oi = e/Qniggs, and $o = 27r/e in the units h = c = 1 
with e and Quiggs being respectively the charge of a test particle and the charge of the 
Higgs particle [H- Experimentally, the situation of an infinitely thin flux tube is realized 
when a flux tube has a radius R which is negligibly small when compared to all other 
length scales in the system. Therefore, both a flux tube with a nonzero and the zero radius 
will be considered. In the case when i? > 0, we shall allow generally for some additional 
interaction inside a flux tube, since such an interaction arises, for example, in the case of 
the magnetic moment coupling. The limit i? — will then depend on the physics inside 
the flux tube. In a rigorous mathematical sense, different physics inside the flux tube will 
be described by different Hamiltonians given as a certain self-adjoint extension of a formal 
differential operator. 

First we shall concentrate on the calculation of the single-particle density of states 
(DOS). The reason is that the DOS is the quantity of basic interest and provides an im- 
portant link between different physical quantities. Knowledge of the DOS determines [via 
the Laplace transform, see Eq. ( |140| ) below] the partition function Z{f3), virial coefficients, 
and in the case of the Dirac equation a relation between effective energy, induced fermion 
number, and the axial anomaly 0, ^. It has been used |Q to calculate the persistent 
current of free electrons induced in the plane by the AB potential The DOS will 
be calculated in two different ways: first, directly through the resolvent and second, by 
only using the scattering properties of the AB potential. After the DOS and differential 
scattering cross sections are calculated various applications are discussed. 

The plan of the paper is as follows. We shall start with the case of the impenetrable 
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flux tube when the flux tube is exterior to the system. In Sec. |^ the basic facts about the 
nonrelativistic AB scattering in this case are summarized. There are neither bound states 
nor zero modes and the DOS is determined solely in terms of the continuous spectrum. 
The single-particle DOS is calculated there directly from the resolvent. One finds that for 
integer a resolvents do differ by a phase factor. However, when the arguments coincide 
they are identical and they do have the same trace and yield a DOS identical to that in 
free space. We shall confirm the anticipation of Comtet, Georgelin, and Ouvry ||^ that 
the change of the DOS is concentrated at the zero energy [see Eq. (^) below]. 

Different self-adjoint extensions in the case of a penetrable flux tube are discussed 
in Sec. ^ Let the coupling constant be negative, e = — |e|. Then different self-adjoint 
extensions arise because, in the channels / = —n and / = —n — 1, where n is the integer 
part of a, there are two linearly independent square integrable at zero solutions of an 
eigenvalue equation. These two channels are also the only two channels where a bound 
state can occur. Bound states are shown to parametrize different self-adjoint extensions 
and the change in the conventional phase shifts. Calculation of the change A; is reviewed. 
Knowledge of the phase shift is then used in Sec. ^to calculate the S matrix and scattering 
cross sections. The S matrix is shown to depend nontrivially on a and not to be a periodic 
function of a. Nevertheless, the S matrix gives rise to differential and transport scattering 
cross sections which are periodic functions of a with period 1. In the presence of a bound 
state, one finds in contrast to the case of the impenetrable flux tube that the differential 
scattering cross section ceases to be symmetric with respect to the substitution — » — 
[see Eq. (0)], where is a scattering angle. The asymmetry is easy to understand 
because (for a > 0) bound states are only formed in the channels for which / < 0. The 
results are then used in Sec. |^ where the validity of the Krein-Friedel formula ||^, ^ 
for the DOS is established for a singular potential. It is shown that the Krein-Friedel 
formula when regularized with the C, function gives the correct DOS. This enables us to 
calculate the DOS for different self-adjoint extensions. The Krein-Friedel formula gives 
the DOS as the sum over phase shifts [see Eq. (|66D] and thereby relates the DOS directly 
to the scattering properties. It is therefore very useful to have its extension to singular 
(especially Coulomb) potentials. One finds that whenever a bound state is present in 
the spectrum it is always accompanied by a resonance [Eq. ([T^D]. Rather surprisingly, 
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Figure 1: Black layer is a superconductor of type II put on top of the heterostructure 
containing the two-dimensional electron gas (shaded region). When this sample is put in 
a homogeneous magnetic field, magnetic field penetrates the superconductor in Abrikosov 
vortices. Therefore, electrons from the heterostructure do not move in the homogeneous 
magnetic field but in the field of a penetrable flux tube. 

the shape of the resonance [given in Eq. (^)] is not of the Breit-Wigner form. Since the 
latter is a direct consequence of analyticity it poses an interesting question on the analytic 
structure of scattering amplitudes for singular potentials. The existence of the resonance 
is quite unexpected. It will influence the transport properties of electrons in the currently 
almost accessible experimental regime |TU|, and the persistent current of free electrons in 
the plane In the limit of a zero energy bound state the resonance goes to zero, too, 
where it merges with the bound state into the continuous spectrum leaving behind the 
phase-shift flip. 

In Sec. P the calculation of the number of bound states for a flux tube of finite radius 
i? 7^ is given. The situation is considered where a short-range potential is placed inside 
the flux tube. The motivation is that vortices that can be realized in real experiments, 
such as in a superconductor of type II, definitely do have a nonzero radius. Another 
realistic realization of the penetrable flux tube is that suggested originally by Rammer 



and Shelankov W^ and later realized experimentally by Bending, Klitzing, and Ploog [ITO 



i. e., to put a type II superconducting gate on top of the heterostructure containing the 
two-dimensional electron gas (2DEG) (see Fig. [^). When a magnetic field is switched on 
the conventional superconductor is penetrated by vortices of flux with a = 1/2. Therefore, 
electrons from the heterostructure do not move in the homogeneous magnetic field but in 
the field of a penetrable flux tube. Electrons can penetrate to their core, in which case the 
potential arises as the result of the magnetic moment coupling. If the gyromagnetic ratio 
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gm is less than 2, < 2, the couphng with magnetic field is not sufficiently strong enough 
to form bound states. If the gyromagnetic ratio gm is exactly equal to 2 (i. e., the magnetic 
moment is not anomalous) zero modes may occur in the spectrum [[T2| , |1^. Their number 
equals, respectively, n — 1 or n, where n is the integer part of a, depending on whether the 
ffux is an integer or not. There is no zero mode for a < 1 ||12|, 0. In the region g^ > 2, 
i.e., exactly where the gyromagnetic ratio of electron {gm = 2.002 32) lies, the coupling 
with magnetic field is enhanced and bound states do occur. In contrast to zero modes |T2| 



one finds that the number of bound states does depend on a regularization. For example, 
in the case of the cylindrical shell regularization |]T1| their number is generally higher than 
for the homogeneous field regularization [see Eq. ( |104D ]. The differences are attributed 



to the different energies of magnetic field inside the fiux tube. In any case, however, the 
number of bound states does not match with the number of zero modes. The question 
about the existence of a bound state in the / = —n — 1 channel is discussed. Although 
generally it is true that only one bound state is present for < a < 1, provided 1 — a is 



sufficiently small, the second bound state does appear (cf. Ref. [0). 



An interpretation of different self-adjoint extensions and the i? — > limit are studied 
in Sec. |^. The existence of a critical coupling is established. In the case of the magnetic 
moment interaction with the interior of the ffux tube, the critical coupling corresponds to 
the case of the normal magnetic moment with the gyromagnetic ratio gm = 2. Provided 
the coupling with the interior of the ffux tube is weaker, bound states do not form at all. 
If the coupling is stronger than the critical value, then, although bound states do exist for 
any 0, they disappear in the limit R ^ 0. Although it might seem surprising at ffrst, 
the fact is a translation of the result of Berezin and Faddeev established more than 



30 years ago that the nontrivial R —y limit of a potential formally given by the Dirac 
delta function 6{r) requires a renormalization of the coupling with the interior of the ffux 
tube (see Ref. for more details). Reference [|18] provides some other examples when 



one encounters the necessity of a renormalization in the quantum mechanics. At ffrst 
sight the renormalization might seem to be artiffcial and only a mathematical obscurity. 
However, one knows that the anomalous magnetic moment always has a form factor that 
does depend on the energy |T^. Therefore, in the realistic situation the renormalization 
is provided by nature itself. 
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The point spectrum [20| at the critical couphng is also considered. One finds that in 



the limit R there are no zero modes in the AB potential for any a. The Aharonov- 



Casher and the index theorems |T2|, |T3|, ^ are corrected in the sense that they do not give, 
respectively, the actual number of zero modes in a given finite-fiux background but rather 
an upper bound. In the presence of a singular field configuration such as the AB potential, 
the square integrability of solutions must also be checked at the position of a singularity of 
the field. It is shown that at such a point the square integrability fails. Another argument 
showing why it happens is to note that the only two channels in the R ^ limit where 
the spectrum can differ from the conventional one are / = —n and / = —n — 1. However, 
for R zero modes [see Eq. (0)] only occur in channels > Z > —n + 1, and hence 
they are never present in those described above. 

In the presence of a bound state the conventional phase shift (|^ acquires a generally 
energy dependent contribution (|37D. In the limit -R ^ bound states are possible in two 
different channels, / = —n and / = —n — 1. However, when R 0, the bound state does 
not occur generally in the / = —n — 1 channel. Therefore, it is natural to expect that the 
phase-shift flip occurs generally only in the / = —n channel. According to this discussion 
the conditions for the occurrence of the phase-shift flip given by Hagen are necessary 
but not sufficient. Moreover, since the origin of the attractive potential inside the fiux 
tube can be arbitrary, our calculation shows [Eqs. (^) and (^)] that the phase-shift fiip 



occurs even in the absence of the spin. Also, in the case of particles with a spin it is 
not the spin but the magnetic moment coupling that is the primary source for the phase- 
shift fiip. A nice interpretation of the phase-shift fiip appears if the bound state energy 
is renormalized to zero. The phase-shift fiip then occurs as the result of merging the 
bound state and the resonance into the continuous spectrum. Although the observation 
of the phase shift fiip is usually attributed to Hagen |T^, it was observed earlier in a 



complementary situation: in the scattering off a general two-dimensional magnetic field 
satisfying the finite-fiux condition in the long- wavelength limit (see Ref. p2[, pp. 437 and 
438). A duality can be observed [see Eq. ( p.ll| )], important from the experimental point 
of view. Imagine two fiux tubes with different radii which are otherwise negligibly small 
when compared to other length scales in the system. Then the scattering properties of 
the fiux tube with a radius Ri at momentum ki are the same as that of radius R2 at 
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momentum k2, provided that kiRi = ^2-^2- 

Starting from the next section apphcations to the total energy of the system (particles 
and a magnetic field), the Hall effect, a persistent current, and the 2nd virial coefficients 
are considered. We shall show that our results are not only of academic but also of 
practical interest thanks to the recent developments in the fabrication of microstructures 
and in mesoscopic physics (see Ref. I^S] for a recent review). The total energy of the 
system and its stability against magnetic field creation are discussed in Sec. ^. In general 
the diamagnetic inequality |]2^ tells us that that the matter is stable. The latter was 
proven under the assumption of minimal coupling. In our calculations we shall allow for 
a speculation that the magnetic moment is independent of momenta. We shall ignore 
the fact that the anomalous magnetic moment has a form factor that vanishes at high 
momenta. Under these hypotheses one can show that in the nonrelativistic case in 2 + 1 
dimensions, a window may exist for the magnetic moment > 2 in which the inequality 
is violated. The reason is the formation of bound states that decouple from the Hilbert 
space by taking away negative energy. Although the dynamics, return fluxes, and the 
form factor must all play an important role in the full (quantum-field theory) discussion 
of the stability, our relation ( |1 14| ) is nevertheless interesting because it gives the stability 
condition in terms of the ratio of the rest and the electromagnetic energies. 

In Sec. ^ we shall examine consequences of the asymmetry of differential scattering 
cross sections. The asymmetry has important consequences as it give rise to the Hall effect. 
The Hall resistivity is then calculated in the dilute vortex limit [see Eq. (|118|) ], i. e., when 
the multiple-scattering contribution is ignored. In Sec. |10| the results are applied to the 
persistent current of free electrons in the plane pierced by a flux tube. Both the spinless 
and the spin one-half cases are dicussed. The above mentioned resemblence between the 
electromagnetic AB potential and the field produced by a cosmic string will enable 
us to conclude that the persistent current might appear in the latter case, too. Once 
again condensed matter physics gives both the motivation and provides a test laboratory 
for a phenomenon that can occur at different length scales in the Universe. Finally, 
in Sec. |Tl|, the 2nd virial coefficient of anyons interacting with a pairwise interaction 
proportional to the Dirac 6 function is calculated. Note that such an interaction appears in 
the nonrelativistic limit of planar field theories p5| , pB| . It will be shown that the 2nd virial 



coefficients are remarkably stable when the interaction is switched-on. The calculation is 
performed directly in the continuum without any use of the customary devices that makes 



the energy levels discrete, such as a finite box or harmonic potential regularization 
0, We shall show that the use of the (^-function regularization reproduces the correct 
answer for the 2nd virial coefficients of noninteracting anyons R, E^, and generalizes 



results of Blum et al. |]28[ for nonrelativistic spin one-half anyons. 

Note that one has the unitary equivalence between a spin 1/2 charged particle in a 
two-dimensional (2D) magnetic field and a spin 1/2 neutral particle with an anomalous 



magnetic moment in a 2D electric field |29]. In our presentation we confine ourselves 



essentially to the nonrelativistic Schrodinger and Pauli equations. The results for the 
Dirac and the Klein-Gordon equations, the induced fermion number, the relation between 
the phase-shift flip and the axial anomaly, and the DOS in the spacetime of a gravitational 
vortex in 2 + 1 dimensions ||30| are discussed elsewhere 0]. For the problems related to the 



gauge transformations that are not discussed here we refer to the review of Ruijsenaars 



2 Impenetrable flux tube and the density of states 

Let us start with the nonrelativistic case and an impenetrable flux tube. We shall consider 
the Pauli Hamiltonian, 

(p - f A)2 
2m 

where fi = fis/s is the magnetic moment operator, s is the spin operator, and s is the 
magnitude of the particle spin. For an electron /Xe = — fl'm|e|/i/4mc = —jiBgm/'^, where hb 
is the Bohr magneton and is the gyromagnetic ratio that characterizes the strength of 
the magnetic moment []I9[. Naively = 2 but one knows that within quantum electro- 



dynamics acquires radiative corrections which depend on the fine structure constant 
c^QED (= 1/137). Within this framework we have [see Eq. (118.4) of Ref. JT^ ] 

eh Q_328^\ _ (3) 



2mc \ 2tt tt 



By separating the variables and assuming e = — |e|, the total Hamiltonian is written 
as a direct sum, H = (BiHi, of channel radial Hamiltonians Hi in the Hilbert space 
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L2[(0,oo),rc?r] § 0, 



Hi = --^--— + — + gm-sJ{r). 4 
dr^ r dr r 



Here = |/ + a|, a is the total flux $ in the units of the flux quantum $o = hc/\e\, and 
S2 = ±1/2 is the projection of the spin on the direction of the flux tube fl], |3l|. The 
Schrodinger equation is recovered upon setting = 0. 

In the case of the impenetrable flux tube the spectra of both the Pauli and the 
Schrodinger equations are identical. Let us first consider the conventional setup where 
wave functions are zero at the position of the flux tube. There are neither zero modes 
nor bound states in this they are incompatible with the boundary conditions. To 

discuss the spectrum, note that for positive (negative) energies the eigenvalue equation in 
the /th channel reduces to the (modified) Bessel equation of the order v = | / + a|, 

Hii^i = k%, (5) 



with k = y/2mE/h. The boundary condition selects only regular solutions at the origin 
and the "spectrum" is given by 

Mr,¥>) = J\i+a\Me''''. (6) 

Generically, to specify the boundary conditions at infinity (or at zero) is necessary only 
for those I for which all solutions of (^ are square integrable at infinity (or at zero): the 
so-called limit circle case (see Ref. p. 152). In general, square integrability takes the 
place of boundary conditions at infinity (or at zero) if one of the solutions of (^ is not 



square integrable at infinity (or at zero): the so-called limit point case (see Ref. p. 
152). From the rigorous point of view, one can speak about an impenetrable flux tube 
only provided that a is not an integer. Then all wave functions ipi are zero at the 
origin. Since Jo(0) = 1, this is not the case for a flux tube with an integer flux and the 
physics may be different from that in free space. In what follows, we shall ignore this 
subtlety, and assume that the / = channel wave function ipi (^) is in the spectrum, and 
even in this case we shall call the flux tube an impenetrable flux tube. 

The AB potential ([I|) is long-ranged and the conventional phase shifts 5;'s [|T|, 

Si = lni\l\-\l + a\), (7) 



are singular, they do not depend on the energy and do not decay to zero for E ^ oo. 
Relation (|^) can be intuitively understood as follows. The AB potential creates "vorticity" 
—a, and positive and negative angular momentum wave functions "go around" the origin, 



respectively, in the anticlockwise and the clockwise directions ||31[| . 

The DOS will be calculated directly from the resolvent (the Green function) (^^(x, y, E+ 
ie) according to the formula 

p^(E) = --ImTrG'«(x,x,E + ie). (8) 

TT 

The integrated density of states Na{E) is then as usual given by 

N^{E)= r p^{E')dE'. (9) 



The eigenfunction expansion for the resolvent in polar coordinates x = (r^;, ({)x) is P3 | 

G«(x,y,E) = -2 / Y: e^'(^^-^^)j|,+„|(A:rJJ|,+„|(A:r,). (10) 

nn Jo Q"^ — , ^ 



By using Eq. (|]), one can check that the resolvent (|T0|) gives the two-dimensional free 
density of states Pq{E) = (m/2nfi^)V for a = n G Z, i. e., when a is an integer, with 
V = J (Pr being the infinite volume. In the latter case, the sum in Eq. ([10| ) can be taken 
exactly by means of Graf's addition theorem (Ref. |3^, relation 9.1.79), 



^ e^'(^-^^) J|,+„l(fcrJ J|,+,|(A;r,) = e-'"(^-^^) Jo(A:|x - y|). (11) 



(=— oo 



By taking the integral in Eq. (|TOD (assuming that = + ze, and < Vy) one finds 



r ~ ^'^ = -^lH^o\q\^ - y|). (12) 



Therefore, 

m 

"2 



G„(x,y,E) = -^-^e-^"(^-^^)/f«(g|x-y|). (13) 



One can also perform the integral 

r .^"^^ ,, Mkr.)Mkry) = ~^MqrM'\qry) (14) 
Jo q"^ + te — 2 

first and then take the remaining sum with the same result (0). Note that Eq. (|T^) is also 
valid for nonintegral u. The Bessel functions are analytic functions and the integrals in 
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Figure 2: Contours of integration in the complex plane. 

Eqs. (|T^) and are taken simply by the residue theorem by using a suitable contour in 
the complex plane (see Fig. ^ and Appendix From the 'formal scattering' point of view 
taking the residue at k = \q\ + ie corresponds to choosing outgoing boundary conditions. 
From the point of view of L'^[{0, oo), rdr] it corresponds to taking the boundary value of 
the resolvent operator on the upper side of the cut at [0, oo) in the complex energy plane 



|35l. Note that 

m 



^+(x,y,i?) = -^—^e^^^^--^y)J^(qr,)Hi'\qry) (15) 



is the resolvent of Hi in the upper-half of the complex plane of q [^. Ji,{qr^) is square 
integrable (in the measure rdr) near the origin, and H^\qry) is square integrable near 
infinity. Obviously, the resolvent is not defined for real positive q since there is no solu- 
tion of the Bessel equation that is square integrable at infinity for such q (see also 
asymptotic formulas in Appendix 0). The point spectrum in a strict sense is empty. 
The spectrum is purely continuous and lies on [0, oo) where the partial resolvent ^^"^ (the 
Fourier transform of the retarded Green function) has a cut The partial resolvent 



Qi in the lower-half oi the complex plane (the Fourier transform of the advanced Green 
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function) is the complex conjugate of P3|] , 



(16) 



Note that the total Green function (^^(x, y,i?) in the presence of an integer flux n ^ 
[see Eq. (^)] differs from Go{'x,-y, E) by a phase factor. However, to calculate the DOS 
one needs the value of G„(x, y,_E') at ipx = and hence the result for the integer flux 
will be the same as at a = n = 0. The limiting value of the resolvent operator on the 



lower side of the cut is the complex conjugate of (|13D , the discontinuity across the cut 



Gn{y^, X, E+) - G'n(x, X, E_) = -i{m/h'^), 



(17) 



and 



- (l/7r)ImGn(x,x,E+) = m/{27ih^) (18) 

which confirms our normalization. 

Whenever a ^ Z Graf's theorem cannot be used. To proceed further with this case we 
use the fact that (p!^ ) has an analytic continuation on the imaginary axis in the complex 
momentum plane 

kdk 



+ k"^ 



Ju(ykrx) Ju(ykry 



J,{iqr^)H^J\iqry) = h{qr^)K,{qry), (19) 



where ly and K^, are modified Bessel functions ||3^. One can obtain this result either by 
performing the integral directly or by performing the analytic continuation in (ITip. To 
sum over / one uses the integral representation of these functions P^ , |5B| . Following the 



steps given in Ref. |^ one can separate the a-dependent contribution 



AG„(x, X, E) = G„(x, X, M) - Go(x, x, M) 



msin(777r) /- r ^^e-Afr.(cosh^^+cosh..) e^^''-'^) 



h' (27r)2 



1 + e' 



(20) 



where M = —i^/2mE /h and rj is the nonintegral part of a, < ?7 < 1. After taking 
the trace over spatial coordinates, using formulas 3.512.1 and 8.334.3 of Ref. ||3^ (see 
Appendix 0), and returning back to the real momentum axis, one finally finds 

1 



TrAG„(x,x,E) = -r/(l-r/) 



2E 



(21) 
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Hence, the change Apa{E) = Pa{E) — po{E) of the DOS induced by the AB potential in 
the whole space is concentrated at zero energy where it is proportional to the 6 function, 

Ap^E) = --IniTr AG4x,x,E+) = -i 77(1 - r^) 5(E). (22) 

TT Z 



This result is similar to the result (|137|) that was obtained in the context of anyonic 



physics (see also Sec. ^ below). The same term has been obtained by Berry as 



the leading term in the semiclassical approximation to the change of the DOS for the AB 
(circular, heart shaped, and Africa shaped) quantum billiards. Thereby, the semiclassical 
approximation turns out to be exact in the R ^ 00 limit, where R is the characteristic 
length of the billiards. 

Let us now take the case where a is negative and of the form a = —n — f] < where 
77, > is a positive integer and t] (the fractional part of a) satisfies < 1] < 1. One can 
readily repeat the same procedure to obtain the change in the DOS. The result is exactly 
the same as in (p2D , 

Ap_,(E) = Ap^iE) = r/(l - 5{E). (23) 

Therefore, Apa{E) as a function of a is a symmetric function under the transformation 
a — >■ —a. Moreover, one can show that Apa{E) as a function of a is only a function of a 
distance from the nearest integer. Indeed, under the substitution 

T]' = l-r] (24) 

one finds that Apa{E) is form invariant, namely 

Ap^iE) = -1 r^(l - r^) S{E) = -\ - i) 5{E). (25) 

3 Penetrable flux tube and self-adjoint extensions 

Physically, self-adjoint extensions in the presence of an Aharonov-Bohm potential can be 
understood as follows. In channels with |/ + q;| > 1, one has a unique choice for the 
resolvent of Hi in L^[(0, cxd), rrfr]. Apart from J\iJ^a\{kr), there does not exist any other 
linearly independent solution of Eq. (^) which is square integrable at the origin. Similarly, 
apart from if|;^^|(A;r), there does not exist any other linearly independent solution of Eq. 
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(|^) which is square integrable at infinity in the upper (lower) complex half- plane of k. In 
this case, the Hamiltonian Hi is said to be in the limit point case at both infinity and 



the origin |3^, and the square integrability takes the place of boundary conditions. An 
ambiguity in defining the resolvent Gf [see Eqs. (|T^) and ([T^)] of Hi arises only in the 
channels with |/ + a| < 1, i. e., with / = —n, —n — 1. There are exactly those values 
of I for which Eq. has two linearly independent (and hence all) square integrable 
at the origin solutions. These solutions can be taken to be J\i+a\ and J-\i+a\- Then, if 
J\i+a\{(irx) is replaced in either (|T5|) or (|T6|) by any linear combination of J\i+a\{(irx) and 
J_|i+a|(gr2.), one obtains a well defined resolvent of Hi. To any particular choice of the 
linear combination corresponds a particular self-adjoint extension of Hi. For these values 
of /, the AB potential is said to be in the limit circle case at zero (see Ref. p. 152), 
and boundary conditions at the origin must be specified to define a resolvent uniquely. 

From the formal mathematical point of view the above discussion is repeated as follows. 
Formal differential operator Hi in [(0, oo), rrfr] is real and hence commutes with the 
involution operator (complex conjugation) in the Hilbert space. Therefore, by the von- 
Neumann theorem (see Ref. p. 143), when Hi is defined as a symmetric operator 
on a dense set in [(0, oo), rrfr] such as C^{0,oo) (the set of infinitely differentiable 
functions that are zero at the origin and decay exponentially at infinity) it always can be 
extented to a self-adjoint operator. According to theorem X.ll of Ref. |3^, except for the 
channels with \l + a\ < 1, the operators Hi are already essentially self-adjoint on C^(0, oo) 
[their deficiency indices are (0,0)]. Only the operators Hi with |/ + a;| < 1, i. e., with 
l = —n, —n — 1, admit a nontrivial one-parametric family of self-adjoint extensions [pi] , P3| 
(their deficiency indices are (1, 1) [^). Therefore, in these two cases the Hamiltonian as 
a self-adjoint operator in the Hilbert space is not defined uniquely. Different self-adjoint 
extensions give different Hamiltonians which correspond to different physics inside the 
flux tube and different boundary conditions at its boundary (see an example in Ref. |p2|| , 
pp. 144 and 145). 

To identify the particular self-adjoint extension one starts with the flnite flux tube of 
radius R and takes the limit i? — > (see Sec. ^. We shall consider the situation with 
bound states of energy Ei = —(fi^ /2m)Kf in the l = —n and l = —n — 1 channels, and with 
n = [a] denoting the nearest integer smaller than or equal to a. This happens when a 
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sufficiently strong attractive short-range potential is placed inside the flux tube. In the 
nonrelativistic case for spin-down fermions this is the case when their magnetic moment 
exceeds the value 2 such as in the case of the electron which has Qm = 2.002 32 0, For 
negative energy, the eigenvalue equation with Hi leads to the modified Bessel equation 
and the wave function is given by 



(26) 



K\i+ct\{.i^if) decreases exponentially for r — > oo and for |/ + a| < 1 it is in i^^[(0, oo), rrfr], 
although it is singular at the origin (see Appendix P). Nevertheless, if we want the bound 
state to be in the Hilbert space the scattering states (|^) in the channels / = —n and 
I = —n — 1 have to be modified. They become a combination of the regular and singular 
Bessel functions, 



ipi{r,ip) = J\i+a\{kr) - AiJ_ii+a\{kr) 



(27) 



This is because Hi has necessarily to be a symmetric operator. This means that the radial 
parts x(^) of any two states ^/'i(r) and ^'2(1') in the Hilbert space have to satisfy 







(2^ 



in the limit r ^ 0, where ] denotes the Wronskian. The condition ( [28| ) is nothing 
but the boundary condition at the origin: in the limit r — > the logarithmic derivative 
rx'(r)/x(r) of any state in the given Hilbert space takes a fixed value This is a 

translation of the mathematical analysis in terms of deficiency indices into more 
physical terms. Obviously, the Wronskian will vanish when the asymptotic forms for 
r ^ 0, up to order 0{r), of any two states are identical. For the general scattering 
solution (^) for < z/ < 1 one has (see Appendix P) 



J,{kr) - AJ^^ikr) 



A 



kr 



1 rri 



AT{l + iy) 



kr 



+ 0{r 



2-u\ 



(29) 



The asymptotic form of ( p6D is then determined by that of K^{z) [see Eq. ( |1 76| ) of 
Appendix |B| 



Kjz) 



Til 



+ 0{z 



2-u\ 



(30) 
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Therefore, the condition (pSf ) determines Ai in the channels l = —n and l = —n — 1 to be 

= {k/^.n?' = {E/\E.n\)\ 

A_^_, = {k/^_^.,f'"'^^ = {E/\E_^_,\Y-\ (31) 

i. e., energy dependent. According to ( P7[ ) and ( pi] ) the bound state energy determines 
the spectrum of Hi in the Hilbert space that is different for different bound state energies. 
Therefore, in physical terms, it is the bound state energy that parametrizes different 
self-adjoint extensions. 

Quite surprisingly, the influence of the bound state on the wave function (|27| ) is most 
pronounced in the limit Eh ^ 0. Then — > oo and the regular wave function changes to 
the singular one, 

^i{r, <f) = J\i+^\{kr) e''^ ^ V^^r, y;) = J-\i+^\{kr) e*''^. (32) 

On the other hand, in the limit of a bound state with infinite bound energy, Ai ^ Q and 

i)i{r,^) = [Jii+a\ikr) - A;J_|,+„|(fcr)] e*''^ ^ 

iJi{r,^) = J\i+a\{kr)e'''^. (33) 



By using the formula (9.2.5) of Ref. [P] one finds that the radial part of the general 
solution (^) behaves for r — oo as follows: 

1-A;e-H'+"1 I ^ ' 

By comparing with its behavior for a = one finds the /th channel matrix to be 
= e^^^^ with 

bi = 5i{E) = - \l + a\) + Ai{E). (35) 

The change A; = Ai{E) of the conventional phase shift (|^ is determined by the equation 

e'^^' = f ..^P (36) 



with the solution given by 



sin(|/ + a|7r) 
^cos(|/ + a|7r) — A^ 



A; = arctan ( ' . (37) 
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Here we have tacitly assumed that Ai as rj ^ 0. For the discussion of this point 
see the end of Sec. |[ Note that when the bound state energy is changed, Ef, — > E'f^, 
and consequently k — > n', then the phase shift 6i{k) in the corresponding channel at 
momentum k equals to S'i{k'), the phase shift in the same channel when the bound state 
energy is E^, provided that 

k' = -k. (38) 

The energy of the bound state breaks the classical scale invariance and sets a scale to 
the problem. Equation (|38| ) then expresses the scaling transformation under which the 
problem is invariant. 

One sees again that the most pronounced influence on the phase shift is in the limit 
Eb t 0. Then A^^ = 0, 

Si = ^7r{\l\ + \l + a\), (39) 

and one has the phase-shift flip when compared to the conventional phase shift (|^) . Using 
the fact that phase shifts are only defined modulo vr one finds 

2 ^ ' 2 2 
5-n-i = -n{n + r]) ~ -na -l^'^^- (40) 

4 The S matrix and scattering cross sections 

Once phase shifts are calculated, the form of scattering amplitudes, the S matrix, and 
cross sections can be discussed. The scattering amplitude /(fc, ip) in two space dimension 
is defined as the coefficient in front of ^^^"^ jr^l"^ in the asymptotic expansions as r — >■ of 
the wave function (see Ref. [HII] for a more rigorous definition). 



where ip' is the direction of the incident plane wave. The scattering amplitude f{k,(f) is 
related to the S matrix, 

iS-l)ik,^) = itk/27Ty/'fik,^), (42) 
where 1 stands here for the unit operator. In the partial-wave expansion one has 

oo 

f{k,ip) = {27nk)-^/^ (^e'^^^W _ l) e^'^, (43) 

l=—oo 
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which expresses the scattering amphtude in terms of the phase shifts 6i{k). The differential 
scattering cross section is given by 



{k,^) = \f{k,^)\'. 



(44) 



In conventional AB scattering, the phase shifts 5;'s are given by Eq. (|^. One then 
finds that 



/ > -n 

/ < - 1, 



(45) 



where, as above, n stands for the integer part, [a], of a. Therefore, the S matrix, s^{(p), 
in the AB potential is given by 

1 



(=— oo 



COS 7ra 



2n 



sm 7r« 



/=— oo 



27r 



— 1 oo 

E -E 

l = -OQ 1 = 



(46) 



where the superscript means that we are dealing with conventional AB scattering. Now, 
the sum in the first term in Eq. (^61) gives the delta function S{ip). Provided the 
second sum in Eq. (|46|) (as well as the first sum) can be taken exactly |^ and one finds 
that the term in parentheses gives 

1 1.1 



-tip 



(47) 



1 - e-'^ 1 - e^^ e'P - 1 
It is obviously singular at = 0. Therefore, care has to be taken to define the sum at 
this point. One has to resist the temptation to extrapolate the result from those at 7^ 
because the result is not continuous at = 0. One can show that the result is given by 



31 



s° ((/)) = cosvra + i 



sm vra 



-mip 



P- 



TT e*^ - 1 • 

where P denotes an analog of the principal value. Formally, 



lim - 

elO 2 



1 



+ 



1 



1 + e 



(48) 



(49) 



The result (^8|) can be verified by taking the inverse Fourier transform and using the 
formula 

e''"^ _ f 71, k>l 
ei<fi - I ~ ] -iT^ k <0. 



dip 



(50) 
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Note that there is nothing particular about the 6 function term in the S matrix in the 
presence of an AB potential. In fact, S{(f) = 6{!f) is the S matrix in free space (in the 
absence of any scatterer), as can be checked directly by substituting a = in Eq. (^). 
The 6 function in Eq. (^) simply means that a fraction cos na of an incident beam passes 
through the flux tube without being scattered. 

In the case of negative a = — |a| the analog of (^5)) is 



e'^^'-l "^^11' ^-^ + ^ (51) 
~ \ e-^-l-l, l<n. 

By repeating the above procedure leading to (^8]), one finds that the S matrix in the 
present case is 

Sin TT rv 1 

s'-H iv) = m cos vr|a| + ^ ^ e^^ " (^2) 

This result for s^^^^{ip) can be verified independently by taking the inverse Fourier trans- 
form and using the formula 

^I/^F-^= -vr, k>0. ^^3) 



A comparison of relations (^) and ( |52| ) shows that, under the transformation a — > —a, 

According to Eq. (^), the scattering amplitude is given by 

/„(A;,<^) = (27rM)V2[s^_l](fc,<^). (55) 



For (p ^ one finds 



f^{k, V) = (l/2vr.fc) V2^^^e-'("+V2).. (56) 

Sm((y9/2) 



Hence, provided 7^ 0, one has by using Eq. (|^) 

/ da°\ , 1 sin^(7ra) , 

The total scattering cross section is infinite for a ^ Z and vanishes for a E Z |^1|. For 
the purpose of the experiment, note that the differential scattering cross section (|57| ) is 
symmetric under (/? —>■ —ip. 
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When bound states are present in the I = —n and in the / = —n — 1 channels, the S 
matrix and the scattering amphtude fa{k,ip) are modified. First one finds that 

sl^sa = sl + Asl (58) 

where 



«=-n-l 



-n 



2i e^^^'+^^'+^'^sinAj. (59) 



l=-n-l 



Here, A; is the change (|37|) of the conventional AB phase shift (^ in the channel /. 
Similarly, one introduces Af^{k,(p) as 

Uk,^) = f,{k,^) + Af,{k,v). (60) 

According to Eq. ( ^Sf ) one has 

— n 

A/°(A;,(^) = (2zM)^/2 J2 e^^^'+^^'+^'^sinA,. (61) 

l=~n-l 



Nevertheless, the property ( [5^ ) of the S matrix still holds and even in the presence of 
bound states one has 

S~\a\{v) = S\a\{-<^)- (62) 

To prove this it is sufficient to show that that a phase shift 6i{E) in the / = —n {I = —n — 1) 
channel for a > is the same as that in the l = n {l = n + l) channel for a < 0. In the case 
of conventional phase shifts one can check for this property directly from Eqs. ( ^5]) and 



(|51|). Since bound state energies and | / + a| remain invariant under the transformation 
a —a and / —I, according to Eq. (0), A/, and hence the phase shift (pS]), remains 
invariant, too. 

Note that if a changes by ±1, the S matrix [see Eqs. ( ^8]) and (^81)] as a function of a 
is not a periodic function of a and changes nontrivially. This is analogous to the fact that, 
as was shown in Sec. ^ the Green's function also depends nontrivially on a. However, 
such as in the latter case when the Green's function gave the single-particle DOS which 
was periodic in a with respect to the substitution a a ±1, the S matrix will be shown 
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to give scattering cross sections that are periodic in a with the same period, too. In the 
absence of bound states this follows immediately from Eq. {^7\j. In the presence of bound 
states, the differential scattering cross section for 7^ is calculated by using Eqs. ( ^41) 
and dl]), 

fc sin(v9/2) ^^^^ - na + (p/2) + sin A_„_i cos (A_„_i + na - (/3/2)|63) 

The periodicity of the differential cross section with respect to the substitution a ^ a±l 
then follows from Eq. (|63D. Note that in the presence of bound states, the differential cross 
section becomes asymmetric with regard to (f ^ —ip (what is equivalent, with regard to 
a —a). The origin of the asymmetry is easy to understand as bound states for a > 
are formed only in channels with / < 0. 

For completeness we also give the result for the so-called transport scattering cross 
section atr, defined by 

{\ — COS (p) — dip (64) 
-TT dip 

[see Ref. Eq. (139.7)]. One finds 

atr = T sm (vra) H — 2^ sm A; 

^ l=-n-l 

— — sin(7ra;) [sin A_„ sin( A_„ — na) — sin A_„_i sin( A_„_i + na)] . (65) 

K 

As is the differential scattering cross section, the transport scattering cross section is also 
periodic with respect to the substitution a a ± 1. 

5 The Krein-Friedel formula and the resonance 

Now, we shall show that to calculate the change of the integrated density of states (IDOS) 
in the whole space, one can use the Krein-Friedel formula |p. The latter gives the contri- 
bution ANa{E) of the scattering states to the change of the IDOS induced by the presence 
of a scatterer, directly as the sum over phase shifts, 

ANo,{E) = - V 6i{E) = (27ri)-i In det S, (66) 
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S being the total on-shell S matrix. As has been shown in Sec. ^, the S matrix (48) is 
singular as a consequence of the singularity of phase shifts (|35| ) which in general do not 
decay when E ^ oo. Therefore, the S matrix (^) or (0) cannot be substituted directly 
into the Krein-Friedel formula (pBf). This can also be seen from the fact that the sum ( pUD 
over phase shifts is not absolutely convergent. We remind the reader that the sum of a 
series that is not absolutely convergent is ambiguous in a conventional sense. If suitably 
rearranged, the (conventional) sum of such a series can take any prescribed value. Thus, 
to deal with such a series, care must be taken. In the present case of the AB potential we 
have found that it is the (^-function regularization which gives the correct answer [ (|22|) ]. 
In the absence of bound states 



oo 

lndetS= J2 '2i5i = t7i I](K|-|/ + a|) 

=— oo /=— oo 

oo oo — n— 1 



a 



1=1 

oo 



l=—n 

oo 



1=1 1=0 1=1 



-i7iri{l-ri), 



(67) 



where and (h are the Riemann and the Hurwitz ( functions (see Appendix 0). Thus, 
by using Eq. (|66D, the change of the DOS is 



(68) 



This result is exactly (p2D. Therefore, despite the AB potential being long-ranged, the 
Krein-Friedel formula can still be used when regularized with the (^-function. In the case 
where a is negative and equals —n — rj, the calculation of the change of the DOS is 
essentially the same. The only change is that now 



I + a, I > n + 1 
—I — a, I < n. 



(69) 



By using the Krein-Friedel formula one again reproduces the result (p3D. 

As has been discussed in Sec. ^, under the influence of bound states the phase shifts 
are changed only in two channels, l = —n and l = —n — 1. Therefore, the contribution of 
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scattering states to the change of the IDOS is then 

A 1 / N 1 / sin(777r) 

ANa{E) = - -7]{1 + - arctan ' ^ ' ^ 



2'' " TT \cosir]n)-{\E^n\/E)^J 

1 / sin(r77r) \ 
— —arctan — r --— , — r , (70) 

where E^n and E^n-i are the binding energies in / = —n and / = —n — 1 channels. By using 



the arguments given in the proof of relation (|62D one finds that the presence of bound 
states does not spoil the property (p3|) of the DOS, and the IDOS is still a symmetric 
function of a, 

AiV_H(E) = AiV|„l(E). (71) 
Note that for < r/ < 1/2 the resonance appears at 

\E I 

Eres = ' > 0. (72) 

[cos(r77r)] 

The phase shift 5_„,(i?) (^) changes by vr in the direction of increasing energy and the 
integrated density of states (^) has a sharp increase by 1. The profile of the resonance 
[the argument of arctan in Eq. (pS])] is given by 

E^tanriTT T 



E^l - E?es E'l - E: 



rj 1 



(73) 



where F is the width of the resonance, 

r = £;;?,,tanr/7r. (74) 



Note that the profile ( [75D is not of the Breit-Wigner form, 

r 



(75) 



E E^-^g 

(see Ref. |^^, §145). For 1/2 < 77 < 1 the resonance is shifted to the l = —n — 1 channel. 

= 1/2 is a special point since resonances occur in both channels at infinity. Therefore, 
the contribution of the arctan terms in ( ffOD does not vanish as 00, but instead 

gives the value —1. Here, we have assumed that bound states are in both the / = — n 
and I = —n — 1 channels. As we shall see in Sec. ^, it is generally easier to form the 
bound state in the / = —n channel than in the / = —n — 1 channel. Therefore, unless r] is 
sufficiently large, only the bound state in / = —n is present and the above considerations 
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-na 
-n/2 



Figure 3: Typical energy dependence of A/(£') in the channel with a bound state. 

have to be appropriately modified. The above discussion also shows that even in the 
generic case, when bound states are present, the DOS still depends only on the distance 



from the nearest integer. Equation ([7^ ) shows that if the energy of the bound state goes 
to zero, Ef, f 0, the resonance also approaches zero energy, Eres I 0. Therefore, in the 
limit Efy ^ the resonance merges with the bound state. As will be discussed in more 
detail in Sec. both the resonance and a bound state disappear from the point spectrum 
and leave behind the phase-shift flip. 

Having calculated phase shifts (|35D for a generic self-adjoint extension, one can calcu- 



late the time delay At associated with the energy derivative of the phase shift [^| 



One sees immediately that if I is different from —n and — n — 1, then the time delay is 
At = 0. The time delay ( [75| ) can be nonzero only in the channels / = —n and / = —n — 1, 
and at the energy corresponding to the resonance the time delay is infinite. 

6 Regular izat ion 

To identify the physics that underlines different self-adjoint extensions we have considered 
the situation when the AB flux tube is regularized by a flux tube of a flnite radius R, and 
the magnetic fleld B inside it satisfles the constraint (pQ]). This situation was discussed 
flrst in Ref. [^. The discussion of this case has a direct relevance for experiment since 
vortices usually realized in experiments are not singular and do have a nonzero radius 
TD| , P^ , |TT|. In order for a bound state to exist, the matching equation for logarithmic 
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derivatives of the exterior (see Appendix 0) and interior (see Appendix |C|) solutions in 
the Ith channel must have a solution. The logarithmic derivative of the exterior solution 
( |26| ) is given by 

Oi{x) = x!^^±^^ < -\l + a\. (77) 

It depends only on the flux a and not on a particular regularization of the interior of the 
flux tube. Here, parameter x is given by 



x = kR= — - R, (78) 
h 

where Eb is the bound state energy. Oi{x) decreases from — |/ + a\ to — oo as x — oo (see 
Appendix |BD. Therefore, 

a(-)< I , (79) 

^ ^ - I a - |/|, I < -n- 1. ^ ^ 

In contrast, the logarithmic derivative of the interior solution depends on both a and 
the particular distribution of the magnetic fleld inside the flux tube. For example, in the 
case of a homogeneous fleld regularization (see, e. g., Ref. and in the absence of the 
magnetic moment coupling [g^ = 0) or any other additional interaction inside the flux 
tube the logarithmic derivative of the interior solution (|186|) is 



^/ N ,,, a M (a + l,b + l,a) 

h{x) = -a+ I +2a- \ \ ' ^ 80 

M [a, 0, a) 



where M{a,b,a) is the Kummer hypergeometric function ||3^, and 

a = — h (x74a) > - 



2 ^ ' ' - 2' 

b = |/| + 1>1. (81) 



For those values of a and b one obtains, by using Eq. (^9_ 

2a M{a + l,b+l,a) 
T M{a,b,a) > 

Hence, 

Ii{x) > -a + \l\ > -\l + a\ (83) 
and it is impossible to get a bound state in this case because the matching equation, 

Oi{x) = Ii{x), (84) 
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cannot be satisfied for any x. Since 



M(0,6,a) = l, (85) 
one finds that Eq. (8^) might have a solution provided a = and a; = 0. Now, our 



purpose will be to break the constraint ( |8lD on a which implies (0). In particular, one 
can show that parameter a can become negative if an attractive potential V{r) is placed 
inside the flux tube, 

nr)|.<.j = -|^-^c(i?), (86) 

and V{r) = otherwise. This amounts to changing x'^/4a in (pl|) to x^/Aa — c/4. If one 
writes c{R) as c{R) = 2[1 + e{R)] one finds 

x^ e(R) 

a = 1 , 

2 4a 2 ' 

b = + (87) 

The attractive potential can be either put in by hand, or, if the Pauli Hamiltonian (^ is 
used, as arising from the magnetic moment coupling of the electrons with the spin opposite 
to the direction of magnetic field B. In the latter case e{r) = const is determined by the 
anomalous part of the magnetic moment, 

e = ig^- 2)12. (88) 

Equation (|87|) shows that the critical value of a = corresponds to the case when = 2, 
i. e., the case when the magnetic moment is not anomalous, and / < 0. In this case the 
matching equation ([8^ ) does have a solution at x = for > / > —n. If one substitutes 



the value of x into either the exterior ( pq) or the interior ( |186| ) solution one would obtain 
a nonsense: "solutions" that do not depend on r (in the case of the exterior solution the 



limit X —> is in fact singular). A more subtle method is needed |T2[ to show that in 
this case zero modes do appear. Since the magnetic field is not singular anymore the 
Aharonov-Casher theorem |T^, |T3| applies. The Aharonov-Casher theorem |T^ tells us 



that in a general finite-flux magnetic field B{r), for which 

/ B(r)dh = $ = const, (89) 
Jn 
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the 2D Pauli Hamiltonian aX gm = 2 has exactly ]«[— 1 zero modes, 

e-<AM^ e-<l'('\xi - 1x2), • • • , e-^^''\xx - ixa)"-^ (90) 

Here, \a \ stands for the nearest integer larger Wioiv or equal to a, and, for a given magnetic 
field -B(r), the function 0(r) is defined by 

0(r) = — / In |r - r'l fi(r') d\' . (91) 
he Jb? 

If a is an integer n, then the number of zero modes is — 1; if not, their number is n = [a] . 



For a < 1 no zero mode is present in the spectrum [0, |T^. The proof of the theorem 
uses the fact that the Pauli Hamiltonian for = 2 can be written as the square of the 
Euclidean two-dimensional massless Dirac operator [O. The result only depends on the 



total flux a and not on a particular distribution of a magnetic field B. The source of the 
attractive potential (RBI) inside the flux tube is not important to the formation of zero 



modes. 

Now, if e > 0, i. e., gm > 2, the parameter a can even become negative. One has 
a < for / < 0, e > 0, and a; G 0, \/2ae^. We shall show that there are at least n + 1 
bound states for any finite R in this case. In other words, if the attractive potential ( ^6]) 
is situated inside the flux tube, the coupling with the interior of the flux tube becomes 
sufficiently strong for the particle to be confined (on the cyclotron orbit) inside it: the 
wave function ( p6D of the bound state decays exponentially outside the fiux tube. In 
what follows we shall confine ourselves to / < 0. The reason is that in order that a be 
negative for / = 1, we must choose gm = 2 + 2e > 6. This is the value of gm that is out 
of experimental interest. Obviously, if Z > 1, then the minimal value of gm that does a 
negative is proportional to / and hence larger than the minimal value of gm for / = 1. 
The number of bound states is given by the number of channels in which the matching 
equation ( P^D can be satisfied with a solution xi > 0. The matching equation ( P^ ) implies 
that the ratio in Eq. (B2]) has to be negative, 



a M{a + l,b+l,a) f 0, I > -n , . 

b M{a,b,a) I < -n - 1. ^ ^ 

This can be satisfied for > / > —n and — 1 < a < only if (see Appendix |C|) 

M{\a\,b,a) <2. (93) 



26 



The latter constraint on the values of a does not impose any physical restrictions since it 



allows for gm € [2,6), i. e., for almost all realistic values of Qm- Since Eq. (|85|) holds, it 
is clear from Eq. (^) that Eq. (|93D can always be satisfied. It is sufficient to look for a 
bound state energy having roughly the form [see Eqs. ( ffSD and 



2t72 

\E,\ = 2aj^{e-t), (94) 

where t (0 < t « e) is some small positive number. Therefore, one concludes that there 
are at least n + 1 bound states at any finite R (cf. Refs. [PUI 0). To satisfy Eq. (P^j) for 



/ < —n — 1 is more difficult because the condition becomes more restrictive. Indeed, it 
is sometimes stated in the literature that for n = and homogeneous regularization only 
one bound state can exist (see, e. g., Refs. |]l5l and |^). We shall show that although 
generally it is true, in some circumstances the {n + 2)th bound state in the / = —n — 1 
channel does appear. It is clear from Eq. (|92|) that the condition to be satisfied gets weaker 
as 1 — ?7 gets smaller. In particular, one always can choose x such that for G (1 — a, 1], 

< (T«l, 

M(|a|,6,a) < 2-7 (95) 

[see Eq. (^5])], where < 7«1 is some small number. Then for these values of x and r/, 
the left-hand side of Eq. (p^) is uniformly bounded, 

a M(a + l,6+l,a) ^ 

I m h \ < -7 < 0. 96 

b M{a,b,a) 

Since a — |/| ^ when rj 1, there exists a o", < a < a, such that Eq. (|92|) is satisfied 
for a given x and rj E (1 — cr, 1) in the l = —n — 1. Therefore, the actual number of bound 
states can in principle be higher than n+1 since the condition ( p2D can be satisfied even for 
some / < —n — 1. An illustrative example is provided by a cylindrical shell regularization 
T^, RUI in which the magnetic field is given by 



i?W = ^5(r-i?). (97) 

As has already been discussed, the exterior solution (26) is the same as in the homogeneous 
field regularization. Only the interior solution changes. We shall denote its logarithmic 



derivative by the superscript 'c'. By using the results of Ref. [^ it can be shown that 
the matching equation takes the form 

Oi{x) = Ii{x) = -a + \l\-ae + ^x^ + 0{x^), (98) 
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where b is as defined by Eq. (|8lD or Eq. (|87|). As in the previous case, the matching 
equation ( pHf ) for Qm > 2 can be satisfied for > / > —n. However, apart from these 
values of /, Eq. ( pH] ) can be satisfied for / below l = —n provided that 

\l\-a<^- (99) 

Note again that when / = —n — 1 one does not generally have a bound state, as is also the 
case of the cylindrical shell regularization. However, if 

l-ne/2 , , 



a bound state does appear in the I = —n — 1 channel. As is seen from Eq. (p9|), when n 
increases, this bound state appears for ever smaller rj. Eventually, for 

n > no =\2/e[ (101) 

the restriction (|99D on rj disappears. Moreover in this case bound states can appear, even 
for a positive /, provided that 



I < (102) 



or, equivalent ly, that 



a > (103) 

However, by taking into account that for the electron — 2 = 0.002 32, the flux has to 
be of order ~ 2000 for this to be the case. To conclude, according to Eqs. (^Uj) and ( |102[ ) 
the number of bound state in the cylindrical shell regularization is 

#b = 1 + n + [a{g^ - 2)/4] + [a{gm + 2)/4 - n] . (104) 

Here [.] denotes the integer part. The number ( [104| ) of bound states in the cylindrical shell 
regularization is generally higher than that in the homogeneous field regularization. One 
can understand the physical origin of this difference in a simple way. In the cylindrical 
shell regularization, the energy Eb of magnetic field is infinite for any R and in this 
sense the magnetic field inside the flux tube is much stronger than, for example, in the 
homogeneous field regularization when 

= f = (lo^) 
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stays finite for any nonzero R. Therefore, in contrast to the number of zero modes given 
by the Aharonov-Casher theorem |13|, the number of bound states does depend not 
only on the total flux a but on a particular distribution of the magnetic field B, and hence 
on the energy of the magnetic field, also. A similar check with the 1/r regularization (see 
Ref. ||3^) allows us to make the hypothesis that their number is less than or equals to 
and that the bound is saturated when one uses the cylindrical shell regularization of the 
AB potential. 

In two dimensions for i? 7^ and A arbitrarily small, the Schrodinger equation always 
has a bound state in the potential XV{r) ( p6D in the absence of the AB potential P3| . 
For R the potential XV{r) is not singular and the wave function of the bound state is 
not singular, either. Now, if the AB potential is put on top of XV{r) the bound state in 
general disappears in its presence. Therefore, the discussion in this section implies that, 
generally, the AB potential has a deconfining effect on the bound state (jTS). Provided 
A > 2 and in the presence of the AB potential the wave function of a bound state ( ^6]) 
is singular and the phase shifts [see Eq. (^)] are changed. Then, in the limit rj ^ 0, 
the singular wave function (^) becomes the regular one and the phase shifts must have 
their conventional values (^. Only the singular bound state wave function necessitates 
the change of phase shifts. This explains why Ai in ( p7D as r/ — > 0, which has already 
been used above in the calculation of A;. 



7 The R ^ limit and the interpretation of self- 
adjoint extensions 

In this section we shall examine the limit R —>■ subject to the condition (^). In the 
case when the flux tube is exterior to the system and particles are not allowed to interact 
with its interior, the R —>■ limit is trivial as there are neither zero modes nor bound 
states. Therefore, in what follows we shall confine ourselves to the case when the flux 
tube is a part of the system and particles do interact with its interior. In the rigorous 
mathemetical sense the limit i? ^ is described by some self-adjoint extension, and we 
shall discuss a correspondence between the i? 7^ case and the limiting case. 

In the limit R the potential V{r) as defined by Eq. goes formally to the 5 
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function, 

a 

V{r)\r<R^-[l + sm 5(r). (106) 

m r 

As has been shown above, until the magnetic moment Qm reaches its critical value Qm = 2, 
nothing is changed with respect to the case of the impenetrable flux tube, and the limit 
-R — s> is again trivial. The limit becomes nontrivial at the critical coupling when 1 
zero modes ( ]90|) exist [1^, [1^ . In this case the symmetry of the spectrum with regard to 



a — s> a ± 1 is lost. Then, as i? we shall show that zero modes ( ]90|) disappear from 
the point spectrum and merge with the continuous spectrum. Indeed, in the limit R —>■ 
the magnetic field B{r) is given by 

B(r) = 2na 6 (r) = -6{r). (107) 
r 

In this case, 0(r) defined by Eq. (|9lD can be calculated exactly. One finds 

0(r)=aln|r|. (108) 

The zero modes (0) are then obviously singular at the origin and they are not elements 
of -^^^[(0, oo), rdr]. As R 0, zero modes ( PO] ) get more and more singular at the position 
of the flux tube and eventually, at the limit, they become nonintegrable and merge with 
the continuous spectrum. In the latter case one has to check for square integrability not 
only at infinity but at the position of singularities of the field, too. It is here where the 
theorems fail. Another argument for the disappearing of zero modes at the R —>■ limit 
is to note that the only two channels at which the spectrum can differ from that of an 
impenetrable flux tube are the channels with I = —n and / = —n — 1. However, for any 
R ^ 0, zero modes never occur in these two channels. Instead, they are in channels with 



> Z > — + 1 [|I2|, [T^. If one suspects that zero modes can appear in some channels 
different from those given in Eq. (0), one can check directly that for whatever /, the 
functions given by Eq. (^0]) are not square integrable either at infinity or at the origin. 
Therefore, in the limit R ^ the symmetry of the spectrum under the substitution 
a — * a ± 1 is again recovered. The best method of illustrating this point is to consider 
a situation when Qm — 2 = 2e > and stays constant as i? ^ 0. Whenever Qm > 2 [or 
Qm = 2 with an attractive potential V{r) = —e/R"^, e > arbitrary small] bound states 
occur in the spectrum. They correspond to solutions > of Eq. ([8^). Note that the 
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solutions are only a function of a. They do not change when K changes. Since x is given 
by Eq. ([78|) one finds that as i? — the bound state energy E^, scales as 

*2 2 

In other words, in addition to the breaking the symmetry a ^ a; ± 1 of the spectrum, in 
the presence of bound states the scale invariance is also broken [^, ^ . Nevertheless, one 
finds that bound states decouple in the R —>■ limit from the Hilbert space L'^[{0, oo), rdr] 
and take away the nonperiodicity of the spectrum (under a —* a±l) that persists for any 
finite R. What is left behind is nothing but the conventional AB problem with the change 
of the density of states given by (|22|) . To show this, note that A/ — > in the general 



scattering solution ([33| ) in the limit k ^ oo. Hence, in the limit the regular state (^) is 



recovered. Another argument is to note that the bound state (p6|) is a function of nr and 
decays exponentially as Kr ^ oo. Therefore, since k ^ oo as i? ^ 0, the wave function 
goes to zero and thereby disappears from the spectrum. The electron is confined inside a 
flux as in a black hole and ceases to communicate with an outside world. 

The above fact might at first be surprising, however it has been demonstrated by 
Berezin and Faddeev |T6| more than 30 years ago that a nontrivial limit -R — > for a 
5-function potential exists only if the coupling constant is renormalized. The latter is 
necessary for a proper mathematical definition of the Schrodinger operator with the 6- 
function potential |T^, . In other words, in order to obtain the bound state in the limit 
-R ^ it is necessary that 

e{R) (110) 

as -R — s> in Eq. (|86D p9| , ^]. The actual energy of the bound state then depends on 
the details of the interaction and the details of renormalization. Obviously, in this case 
the symmetry of the spectrum under a — > a ± 1 is broken. Since the spectrum and phase 
shifts do not change until the attractive potential inside the flux tube is renormalized to 
its critical strength and a bound state is formed, this generalizes the result of Ref. |^ 



that the AB scattering in the case of open boundary conditions at the flux tube boundary 
coincides with the AB scattering with Dirichlet boundary conditions. In other words, 
provided the attractive potential inside the flux tube is not renormalized to its critical 
value and remains either weaker or stronger than the critical potential, then the fiux tube 
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remains impenetrable in the limit R 0. 

If a bound state is present, the phase shift in a corresponding channel acquires an 
energy dependent term (0). The latter, in the limit of zero bound energy, gives rise to 
the phase-shift flip. Therefore, it is natural to assume that the change (|37D of the phase 
shift or the phase-shift flip will occur in the limit R —>■ only in these channels where 
the bound state occurs for i? 7^ 0. When the phase-shift flip takes place, the symmetry 
a a ± 1 is again broken. Our calculation [see relations ( pSj) and (|39D] shows clearly 
that the phase-shift flip |T^ is not connected to the spin but may occur in its absence as 



well. 

As has been mentioned, solutions of the matching equations are only functions of the 
flux, a and x{= kR). On the other hand, scattering solutions are functions of the flux 
a, the ratio k/n [see Eq. (|3T|)], and kr. From the experimental point of view it is useful 
to remark the following duality: the physics at a given radius i?i of the flux tube and at 
momentum ki is identical to that at R2 and k2, provided 

k2 = ——- (111) 

When Eq. ( |1 1 1[ ) holds, then the relative combination of the regular and the singular Bessel 
functions in Eq. (P?]), and hence the phase shift, are the same. Therefore, provided one 
has only a vortex of a finite radius -Ri at disposal, one can examine the physics of almost 
singular vortices with a radius R2 <sC Ri by performing experiments at very small momenta 
k, i. e., such that Ri<<cl/k. It is in the latter situation where the phase-shift flip has 



been established [^. Moreover, it is easier to realize experimentally. 



8 Energy calculations 

In this section we shall assume that Qm is a fixed constant and that no renormalization of 
gm occurs. The energy of the system consisting of particles and field will be evaluated for 
a sequence of flux tubes of decreasing radii subject to the constraint that the total flux $ 
( |89| ) is the same in each. Therefore, dynamical phenomena such as the induction of the 
electric field and return fluxes will be ignored. One reason for this rough approximation 
is that we do not know better. 

Let us now discuss the cases where gm is, respectively, less than, equal to, or greater 
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than 2. It has been aheady shown that up to = 2 no bound state is present in the 
spectrum and the change of the density of the scattering states is still given by Eq. (0). 
Zero modes which occur for g^, = '2 at R ^ are regular at the origin and do not change 
phase shifts as i? — 0. Because the energy of the magnetic field (|105|) tends to infinity 
as i? ^ the system consisting of particles and field is definitely stable with respect to 
spontaneous creation of the AB field. 

When Qm > '2, then bound states occur in the spectrum. Their energy is given by 
Eq. ( [I09| ) and scales to —00 when i? ^ 0, in the same way (as 1 /R^) that the magnetic 
field energy ( |105D does to +00. As has been shown in Sec. ^ provided e{R) is not 
renormalized in the limit -R ^ 0, the bound states decouple in the limit from the Hilbert 
space -^^^[(0, 00), rdr] [Ai — > in Eq. ( p7| in the limit) and take away the nonperiodicity oi 
the spectrum with regard to a — a ± 1 that persists for any finite R. What is left behind 
is nothing but the conventional AB problem with a change of the density of states (|2^) . 

In what follows the homogeneous field regularization will be used. Note that the 
homogeneous magnetic field optimizes the energy functional 

E = [ B'^{r)dh (112) 
Jn 

subject to the constraint (^9|) . Bound state solutions xi for the homogeneous field regu- 
larization determine the function X{a,gm), 

X{a, Qm) = (47ra2)-i ^ x^{a, g^) > 0. (113) 

I 

By comparing the coefficients in front of 1/R^ in Eqs. ( |105| ) and (|109|) , one finds that 
whenever the ratio of the rest energy to the electromagnetic energy is less than X{a, gm), 

2 

TTIC 

^<Xia,gm), (114) 

the total energy of field and matter together goes to —00 as i? 0. Therefore, in the static 
approximation, without the account of the electric field energy, the energy of the system 
decreases with decreasing R. This does not show the instability against the spontaneous 
creation of a magnetic field yet, since the full treatment has to take the dynamics and the 
magnetic moment form factors into account. Nevertheless, the discussion shows that the 
case of gm > 2 is different with respect to gm 2. In three space dimensions, the function 
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X{a,gm) is replaced by X{a, gm)/L, with L the length of the flux string, and one has to 
discuss the density of states for a long flux ring . 

Note, in passing, that in the relativistic quantum-mechanical treatment one flnds 
that the system is stable. The reason is as follows: for spin-up electrons the magnetic 
moment coupling introduces a repulsive interaction and, hence, there is no natural way 
to obtain a spectrum with bound states as there was for spin-down electrons in which 
case this interaction is attractive. To obtain a bound state in the spectrum, an attractive 
potential 2V{r) has to be put inside the flux tube by hand. This is the principal reason 
why in the case of the Dirac electron (a "pair" of the spin-up and spin-down Schrodinger 
electrons) the magnetic moment coupling cannot produce a bound state in the spectrum 
no matter how large or small the magnetic moment is [Q. However, there is one exception 
and this occurs when the gauge fleld is the Chern-Simons fleld [0. Indeed, it has recently 



been discussed by Hosotani |^2[ that in the full-fledged quantum-fleld theory model with 
the Chern-Simons gauge fleld a magnetic fleld can be spontaneously generated. The 
Chern-Simons fleld is somewhat pathological with respect to the discussion in this section, 
for in this case the density of matter acts as the source of the AB gauge fleld, and all 
particles carry a flux In this respect one has "spontaneous magnetic fleld generation" 
whenever the particle density is different from zero. The discussion that led us to the 
stability condition (|1 14| ) remains reasonable even in this case. Since the energy of the 
Chern-Simons fleld is zero, there is nothing to impede the formation of a magnetic fleld. 
Regarding massless charged particles, note that the result of Gribov concerning an 
instability of massless charged particles shows that the ratio of the rest energy to the 
electromagnetic energy is an important parameter in fleld theory, and a condition similar 
to ( |114D must hold. The massless charged particles were claimed by Gribov not to exist 
in nature, since they are completely screened locally in the process of their formation. 



9 The Hall effect in the dilute vortex limit 

As has been discussed in Sec. |, the differential scattering cross section ( |63D for a generic 
self-adjoint extension is asymmetric with regard to —>■ —ip. One can show that the 
asymmetry can have important experimental consequences: in contrast to the conventional 
symmetric differential cross section (p7D, the asymmetric differential cross section ( |55D 
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Figure 4: The differential scattering cross section da{k, ip) gives the fraction of an incident 
current J which is scattered to the angle ip. In the case of an asymmetric differential 
scattering cross section there is generally a surplus of particles flowing in one of the 
transverse directions. 



can give rise to the Hall effect. Indeed, if the incident wave function is normalized to 
unit current density, the differential scattering cross section, da{k,ip), is nothing but the 
transition probability between the incident scattering state and the scattering state that 
propagates in a direction ip with respect to the incident wave In other words, the 
differential scattering cross section gives the fraction of particles from the incident beam 
that are scattered off to the angle if. Now, let us consider electrons propagating with 
the Fermi momentum kp in a sample in a direction singled out by an applied electric 
field. If vortices are randomly distributed throughout the sample the electrons will be 
scattered. In what follows, the dilute vortex limit will be considered, in which the multiple- 
scattering contributions are neglected. Results concerning the Hall effect are then obtained 
by summing over the single-vortex contributions. The asymmetric differential scattering 
cross section ( |63D of an electron from a vortex means that, generally, there is a net surplus 
of the electrons propagating in one of the transverse directions, i. e., either to the right or 
to the left (see Fig. ^). The quantity that measures the fraction of the electrons moving 
in a transverse direction is sin(y9) dalkp, Therefore, if the density of vortices is ra^, the 
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Hall current, in the dilute vortex limit, is proportional to 

1 da 

Uy / dip siuLf ——{kF, (115) 

J-TT dip 

By inverting the conductivity tensor one finds that the Hall resistivity, pxy, is 

kp dip da ^ 



Pxy = pH— 7^sin(^— (/Ci.,(^), (116) 

a J-vr/vr dip 

which was obtained by Nielsen and Hedegaard [|^. In the latter case the result (|116|) 
was obtained by solving (in the dilute vortex limit) Boltzmann's equation, which relates 
the scattering and the transport properties. Here, p^/ = Bc/nf.e is the Hall resistivity in 
a uniform magnetic field B, rie is the density of electrons, and e is the electronic charge. 
The uniform magnetic field B is obtained by averaging over the field produced by vortices 
with the density n^, 

B = n^a^o- (117) 

Fortunately, because of simp factor in Eq. ( |116| ), the differential scattering cross section 
is only needed for 7^ to determine the Hall resistivity. When the differential scattering 
cross section (|63| ) is inserted in Eq. (|116|) one finds that only the last term contributes 



and 



4:11 he? 

Pxy = — ^— ^ sin(7ra) [sin A_„ cos(A_„ — -no) + sin A_„_i cos(A_„_i + Tra)] . (118) 



Equation ( |118|) shows that one needs more than the asymmetry of the differential scat- 
tering cross section for the Hall resistivity to be different from zero. In fact, the Hall 
resistivity, p^y, vanishes whenever (modulo vr) 

A_„ = -A_„_i. (119) 

However, as has been discussed in Sees. H and |^, in any realistic situation relation ( |119| ) is 
not generic and the Hall effect will appear. For example, in the case where the phase-shift 
flip occurs only in the / = —n channel, i. e., A_„ = rjir and A_„_i = 0, one flnds 

An hc^ 

Pxy = — sin(7rQ;) sin(r77r) cos(7rn). (120) 

Tig € 

In this case it is easy to check that if the vorticity a increases, p^y does not change its sign. 
Our result (|118|) shows that the Hall resistivity is proportional to the density of vortices 
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and depends on their vorticity via trigonometrical functions. As a self-consistency check, 
the Hall resistivity ( |118|) vanishes whenever a is an integer. In the case of vortices in a type 



II superconductor [jTO| , if the magnetic field increases, the vorticity of each vortex remains 
constant, and only their density changes linearly with the applied field. Therefore, 
the dependence of the Hall resistivity, p^y^ on the magnetic field is linear in the dilute 



vortex limit [|Tl|]. Obviously, if the magnetic field is sufficiently large, the dilute vortex 
limit ceases to be valid and deviations appear [|1^]. Nevertheless, measurements of the 
Hall effect on a single isolated vortex is now almost experimentally possible |T^, |T3[, and 



the above results for the Hall resistivity can be tested. Note that, provided the resonance 
(JT^) is close to the Fermi energy Ep^ an interesting effect may apear because the Hall 
resistivity becomes very sensitive to the changes of Ep. 

10 Persistent current of free electrons in the plane 
pierced by a flux tube 

The persistent current in a finite (ring) geometry was first discussed in Ref. ||5^. It is 
reminiscent of the edge currents (see discussion in Ref. on their existence) that arise 
in the presence of a magnetic field in the absence of the electric force. In the case of a 
bounded system, the energy levels are discrete and the persistent current induced by a 
flux tube with flux carried by the jth eigenstate, is 

V H — 



= -^^^ ■ (121) 



In the case of an unbounded system the spectrum will have both a discrete and a continu- 
ous part. The persistent current carried by an isolated eigenstate (from a point spectrum) 
is still given by formula ( |121| ). The contribution of scattering states to a persistent current 
is then determined by the formula 

dI{E, a) = {2my^d^ [In det S(^, $)] dE, (122) 

derived by Akkermans et al. 0. Here S(i?,<l') is the on-shell scattering matrix in the 
presence of a flux tube, and dI{E, $) is the differential contribution to the persistent 
current at energy E. The persistent current was defined with respect to a point. It was 
given by the total current through a line that extends from that point to infinity, in the 
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absence of currents through the external leads. Now, by the Krein-Friedel formula (|6^), 
Indet S{E, $) is directly related to the change ANa{E) of the IDOS, and hence (by using 
$ = a$o) "we have 

dI{E, a) = -^da [ANa{E)] dE. (123) 

As has been shown in Sec. |], ANa{E) is a symmetric function of a. Therefore, a persistent 
current is an antisymmetric function of a (see Figs. |^ and ^). One can show that the 
formula ( |123|) reduces to (|121|) in the case of the discrete spectrum, and in fact, the formula 



is valid for both continuous and discrete parts of the spectrum. In the latter case, the 
DOS is formally given by p{E,^) = — Ej{^)], where the summation is over all 

discrete levels. Hence 

*§i)^-i:.'iB-^,(t)i§. (124) 

j 

Therefore, 

J UM/ j,Ej<E 

Now, by substituting the result into Eq. (|123|) and after integrating up to the Fermi energy 
Ep, one recovers the sum over all contributions of single levels, as given by the formula 



12 1|) , below the Fermi energy Ep. 

For spinless fermions neither bound states nor a phase-shift flip occur and the change 
in the DOS is given by (^2]) . Therefore, when all states below the Fermi energy Ep are 
occupied, one finds that the persistent current of spinless fermions around the origin, 
which is pierced by a flux tube, is 

r eEp 



he 



iv-l/2). (126) 



The current depends linearly on t] (cf. Ref. where it is a constant) as it does in small 
one-dimensional metal rings In contrast to Ref. 0, the current is antisymmetric not 
only about the values a = n, where n is an integer, but also about values a = n + 1/2 
where it vanishes (see Fig. |^). The latter values of a are such as the former the values of 
a where time invariance is preserved. 

In the case of spin one-half fermions, the contribution of spin-up fermions to the 
persistent current is still given by formula (|1 26|) . The contribution of spin-down fermions 
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Figure 5: Characteristic dependence of a persistent current of free spinless fermions in the 
plane on the flux a (arbitrary units in the y direction). 

depends on their magnetic moment gm- At the critical value of gm = 2, either the phase- 
shift flip or a bound state can occur. The most general expression for the contribution of 
scattering states to ANa{E) for E > 0, which includes the situation where bound states 
or a phase-sfhift flip are present, is given by Eq. ([7^). The persistent current in this case 
is obtained by substituting the result ([7^) for ANa{E) directly in Eq. ( |122| ). Here, one 
must not forget that the bound state energies also depend on flux (JT^. In the case that 
the phase-shift flip occurs in the / = —n channel, the contribution of spin-down fermions 
to the persistent current is given by the formula 

r eEp 



he 



(127) 



[see Fig. P(a)]. Therefore, the total persistent current of spin one-half fermions in the 
plane when a phase-shift flip occurs is 



he 



:i28) 



[see Fig. P(b)]. This is another important difference from Ref. where they obtained 
a result that the total current is zero. The result is consistent with general requirements 
[01 of periodicity with regard to a ^ a ± 1, and antisymmetry with respect to a — —a. 
Provided that the phase-shift flip occurs both in the / = —n and / = —n — 1 channels, 
the contribution of spin-down fermions to the persistent current is the same as that of 
spin-up fermions. However, as has been extensively discussed in Sees. ^ and |^, the latter 
case of the phase-shift flip in two channels is less probable than that of the phase-shift 
flip in a single channel. 

Observations of the persistent current may finally reveal the resonance ( |72D in the AB 
scattering, since near it the current becomes very sensitive to the change of the Fermi 
energy Ep and of the flux. The similarity between the scattering in the presence of the 
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a) 




Figure 6: (a) The contribution of spin-up electrons to the persistent current as a function 
of the flux a, provided that the phase-shift flip occurs in a single channel, (b) Persistent 
current of spin one-half fermions as a function of the flux a. 
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AB potential and in the field of a cosmic string naturally suggests that a similar current 
should occur in the field of a cosmic string, too. 



11 The 2nd virial coefficient of nonrelativistic anyons 

Anyons are usually represented as either bosons or fermions threaded by the flux tube 
with the flux $. Noninteracting anyons are described by the Hamiltonian 



H 



2m 



(129) 



where is the number of anyons, and Ajj = A(rj — rj) is nothing but the AB potential 
d^) centered at the position of the jth particle. 
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(130) 



written here in a slightly different form with z the unit vector in the direction of the flux. 
In what follows we shall consider anyons in the presence of a pairwise interaction U (r — r'). 



f/(r - r') 



-A a 

m 



5(2) (i 



-A- 



m r 



(131) 



By transforming as usual to the center of mass and relative coordinates (r, ip), and leaving 
aside the free motion of the center of mass, the relative Hamiltonian takes the form 



rel 



m 



dl + -dr - \i-id^ + a 



+ U{r 



(132) 



where, as above, a = $/$o- The form of the relative Hamiltonian Hrei corresponds |^ 



to that used in Refs. |27] and [2^. One neglects the electrostatic forces between anyons 
by assuming the limit e with a fixed. The relative wave function is parametrized as 
exp{iL(f) f{r), where L is standard, namely L = 21 for bosons and 2/ + 1 for fermions. 
After parametrization, the relative Hamiltonian Hrei takes a form similar to that of Hi (|^), 
provided the substitutions X = Qm and Sz = —1 are made, and the reduced mass rh = m/2 
is used. Due to the parametrization of the relative wave function, the parameter v is now 



V = 2l + a 
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(133) 



provided one starts from the bosonic end, or 



z/ = 2/ + 1 + a 



(134) 



when one starts from the fermionic end 

The equation of state of a real gas expanded in powers of the density, /i = N/V , is 



PV 



I + a2/i + as/i + 



(135) 



where the aj stand for the virial coefficients (see, for example, Ref. |^). Here, P is 
the pressure, V is the volume, and (3 = 1/kT. The calculation of a2{T) only requires a 
knowledge of two-body interaction. We shall show that the results of preceding sections 
(with a slight modification) can be directly applied to the calculation of the 2nd virial 
coefficient a2{T) of the gas of anyons. Let us first consider the noninteracting case A = 0. 
By using the Krein-Friedel formula one can calculate the change of the DOS in a way 
similar to that used in Sec. |[ In the case of bosons one finds 



lndetS= J2 225i = «7r ^ (|2/| - |2/ + a|) 

l=—oo l=—oo 
" oo oo oo 

.1=1 1=0 1=1 

= 2m [2Cr(s) - CHis,r]a) - Cnis, 1 - r]a)]\. 
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and hence 

Ap^(E) = -r^,(l-r/J5(E). (137) 

Now, in the case of anyons, rja is the fractional part of a/2. Note that formally, if the sum 
in Eq. ( p?] ) is restricted to even / only, the result ( |137] ) is twice that in the AB potential 
with unrestricted I. In the case of fermions (here rja is supposed to be within the range 

\Va\ < 1/2), 



IndetS = 2^5; = ^ri ^ (|2/ + 1| - |2/ + l + a|) = 

l=—oo l=—oo 

oo oo oo oo 

Y^il + 1/2)-^ + $:(/ - 1/2)- - ^(/ + va + 1/2)-^ - -Va- 1/2)- 
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2m 7]^ 
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and hence 

ApiiE)=r^l6iE). (139) 
Now, the two-body interaction partition function Zint{(3) can be calculated from 

Z':^M = le-^^ApyiE)dE. (140) 

Note that Zint{l3) vanishes for 77a = when interactions (including the AB interaction) are 
switched off. The integration here runs over the whole spectrum. However, since there 
are no bound states, Apa{E) is zero for E < and the integral reduces to the Laplace 
transform. By inserting Eqs. ( |137D and ( |139| ) into Eq. ( |14(j| ) one finds that the partition 
functions do not depend on the temperature T [0, 



int 



-Vail -Va), 



7f 



(141) 



a2{T) can be directly expressed in terms of the two-body partition function and has the 
form 



4(T) = 
2^ ' 2m 



1 + 



int 



in the case of bosons, and 



4{T) 



2m 



(142) 



(143) 



in the case of fermions [0. The prefactor in the last two expressions can be written as 
A|/4, where At = ^j27iP/m is the thermal length. Now Eqs. (|l4l|), ([l42D , and ([l43D 
imply that the 2nd virial coefficients are periodic with respect to rj^ —>■ rja ± 1, i. e., with 



respect to a ^ a ± 2 ||2^ . In case of bosons one finds for < a < 2 that 



4(T) 



A 



1 - 4a 2a' 



A 



2(a 



4 V '74 
Similarly, in the case of fermions one obtains for — 1 < a < 1 that 

4(T) = ^(l-2a^ 



A 



l-4(a + l) + 2(a + l) 



(144) 



(145) 



The result for other ranges of a is obtained by using the periodicity. The 2nd virial 
coefficients o/^iT) and a^iT) are written in two equivalent forms. The second form clearly 
shows that if a is raised to a -|- 1, then 03 (T) — > aliT). Similarly, if a is lowered by 1, 
then 4(T) a^(T). 



43 



So far, we have only reproduced the known results for the 2nd virial coefficients of 
noninteracting anyons [0, In the presence of an interaction (as may be the case for 
anyon-antianyon interaction) the virial coefficients will change. The results for the DOS 
in the AB potential then enable us to calculate them in the special case of the potential 
U{r — r'). The analysis is similar to that of Sec. ^ Therefore, provided the parameter 
A 7^ 2 in Eq. ( |131| ), no bound state is present in the spectrum and the 2nd virial coefficients 
a^'^T) are still given, respectively, by Eq. ( p,44| ) or Eq. ( |145| ). They exhibit quite a lot 
of rigidity with respect to the interaction and they start to change only at the critical 
coupling Ac = 2 when the phase-shift flip or a bound state may occur (at an energy that 
depends on the details of the limit when the radius of the flux tube shrinks to zero). The 
only change with regard to the AB scattering discussed in previous sections is that the 
parameter u changes here by 2 and hence a bound state or the phase-shift flip can only 
occur in a single channel. Then the change of the IDOS ANa{E) is given either by 

AN^iE) = M{va) + - arctan f . f^^.^^^^ ) (146) 

TT \ cos(?77r) - {\E_n\lE)'^ 



or 



where either N{rja) = —Vai^ — Va) or Mirja) = n^ depending on n and whether one starts 
from the bosonic or the fermionic end. To calculate Zint{l3) in this case one integrates 
by parts the general formula ( [L4(J| ) and rewrites it in terms of the change ANa{E) of the 
IDOS [0. Now a bound state is present and one obtains in general 



POO 

ZintiP) = E e"^""" + f3J^ e-^^'AN^iE) dE, 



:i48) 



b 

where the sum here in principle runs over all bound states. Note that at the critical 
coupling, the partition function depends on T. Eventually, the 2nd virial coefficients are 
obtained by inserting the result in either (|144|) or ( p.45| ). In Table I the results are presented 
for the partition function, and the 2nd virial coefficient for bosons and fermions in the 
case when the phase-shift flip occurs for n = [a] irrespective of whether odd {n = no) 
and even (n = ng). The results are presented in terms of 77. Provided n is even, rj is 
added to (subtracted from) ZmtiP) in the case of bosons (fermions). The reason is that 
in the former case the condition < |z^| < 1 is satisfled for the only / given by 21 = —n 
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{21 = —n — 2), which imphes that the parameter v takes the value 7] [—(1 — rj)]. When n 
is odd, the situation is reversed and rj is subtracted from (added to) ZintiP) in the case of 
bosons (fermions). The calculation can be performed straightforwardly, but care has to 
be taken with regard to the range of rja- < rja < 1 in the case of bosons and |?7a| < 1/2 
in the case of fermions. 



TABLE I. The partition function Zint{f3) and the 2nd virial coefficient for bosons and 
fermions in the presence of the phase-shift flip for n = [a] odd and even. 



n 




4{T) 




a{{T) 


even 


i(2r/ + r/2) 


f [l-2(r/ + in 


i(-4r/ + r/2) 


f (1 + 877 - 27^2) 


odd 


-i(l + 4r7-r/2) 


f (1 + 8r/ - 27^2) 




f [l-2(r/ + in 



The results presented in Table I, apart from the periodicity of the 2nd virial coefficients, 
demonstrate clearly that the 2nd virial coefficients interpolate between the fermionic and 
bosonic case as a changes by unity. Namely, the 2nd virial coefficient in the UoB case 
(bosonic end with n = rio odd) is the same as that in the rieF case (fermionic end with 
n = rtf. even), and similarly the result for the 2nd virial coefficient in the rtoF case is the 
same as that in the UeB case. The reason is that in the respective cases, the partition 
functions Z^^^ receive the same contribution. This is not restricted to the case when the 



contribution is given by the phase-shift flip only, but holds in the general case [( |147|) and 
(|148| )] when the bound state is present. Then, according to Eqs. ( p.42| ) and (|143|) , the 
changes in the 2nd virial coefficients are identical. Therefore, the interpolation between 
different statistics as a — a ± 1 holds both with and without the short-range interaction 



( |131| ). There is, however, one distinguishing feature of the 2nd virial coefficients at the 
critical coupling observed in Ref. [^, namely, their discontinuity. Since a discontinuity 
in the virial coefficient implies a discontinuity in the free energy, to the same order of 
approximation this points to some kind of phase transition It has already been 

mentioned that at the critical coupling Ac = 2, the scale invariance is generally broken 
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and one expects a different physics in this case. 



Our results cover naturally the case of the spin one-half anyons |2^, too. In the latter 
if we ignore the anomalous part of the magnetic moment, but add a Zeeman interaction, 
an anyon-anyon potential U{r — r') arises which has the form 

U{r - r') = -(si + S2)^- 5(|r - r'|), (149) 
m r — r 



where Si = ±1 and S2 = ±1 are spin projections on the direction of the flux tube |^ 
Provided that at least one of the spins does not have an orientation opposite to the 
orientation of the magnetic field, then, depending on the statistics, the result for the 2nd 



virial coefficient is given either by Eq. ( |144| ) or Eq. ( |145| ), in accord with and with 



previous results ^7j: the interaction does not have its critical value and one should 
recover the noninteracting case. If both spins are opposite to the direction of the magnetic 
field, the critical potential with Ac = 2 arises and, provided the phase-shift fiip takes the 
place, the results for the 2nd virial coefficient of anyons can be read off from Table I. 
Our result in the HoF case agrees with that of Blum et al. [^. In the UeF case they 
did not take into account that the phase-shift flip may occur as they assumed that in the 
AB scattering the phase-shift flip cannot occur in the l = —n — 1 channel. Therefore, not 
surprisingly, in the absence of the phase-shift flip in Ref. p8[ they obtained the same result 
(|145| ) for the 2nd virial coefficient in the rieF case as in the absence of the interaction 
f/(r — r'). As has been discussed in Sec. ^ it is more difficult to create a bound state in 
the / = —n — 1 channel, and the phase-shift flip in this channel may possibly not occur. If 
this happens, then our result for the 2nd virial coefficient coincides with that of Ref. |^ 



For the case of an unpolarized system, the complete a^iT) is obtained by averaging 
over four spin states. Provided that there is no phase-shift flip in the n^F case, a^iT) is 
the same in all four cases and is given by Eq. ( |145| ). In the presence of the phase-shift flip, 
one of the values 02 (T) is then given by Table I. In the rioF case, three values of a^iT) 
are the same and given by Eq. ( |144| ), because of the interpolation between fermionic and 
bosonic statistics. The remaining value of 02 (T) under the presence of the phase-shift flip 
is again taken from Table I. Finally, one obtains 

^2 C (-'^ ~ 2?7^), UeF case without the phase-shift flip 

a^(T) = ^{ {I + 27]- 2r/2), UeF case with the phase-shift flip (150) 
[—l + 27] — 2rf)^ UoF case. 
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It is worth noting that anyons with a short-range pairwise attractive interaction occur 
in the nonrelativistic hmit, taken up to the terms of order f^/c^, of the topologically 
massive planar electrodynamics |51| defined by the action 



S 



1 

Z2 



matter 



0/3 



a^'o^aA dt(fv 



(151) 



with a nonminimal coupling [25, 26 



J matter 
N 



E 

a=l 



oO 

'-'matter 



dtd r. 



(152) 



Here, |-R| represents the screening length of the electromagnetic interaction, L is the 
arbitrary scale parameter, stands for the ath particle's charge in units of L, J^'^^ denotes 
the standard one-particle current normalized to the unit charge, and is the number of 
the pointlike particles. In the model the photon is massive with a mass M = h/c\R\. Note 
that F^^ is actually a magnetic field, and the nonminimal coupling term that is peculiar to 



2 + 1 dimensions p5|, ^ is in this instance the familiar Pauli magnetic moment coupling 
that exists here even in the absence of a spin. It has been shown by Stern ||25| that, 
provided 

ga = L'R^a, (153) 
§(jdx^Ai becomes a topological invariant. 



doc 



c 



gn. 



(154) 



The screening length remains nonzero, and the electric charge and magnetic moment 
balance each other in such a way that the radiation is absent. At the tree level, the theory 
at the critical coupling ( |153| ) can be exactly reexpressed as a simple effective action at a 



distance without the usual complications associated with the retardation, etc. [^. The 
situation is analogous to that in pure Chern-Simons (CS) theory, 

Scs = ^l t^p-yA^d^AUt dh (155) 

which is recovered in the limit — oo, i? ^ 0, keeping fixed the parameter 

= L'^R. (156) 
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In the latter case the absence of radiation is simply because the photon gets infinitely 
heavy. At the tree level, the theory at the critical coupling ( |153| ) can be exactly reexpressed 
as a simple effective action at a distance without the usual complications associated with 
the retardation, etc. p5 . 



12 Discussion of the results and open questions 

The nonrelativistic scattering in the AB potential has been analyzed. The DOS and scat- 
tering cross sections have been calculated and various applications have been discussed. 
Despite the fact that the single-particle Green function is not a periodic function of a, it 
gives rise to the DOS which is a periodic and symmetric function of the fiux and it depends 
only on the distance from the nearest integer. It has been shown that the Krein-Friedel 
formula |^ is not restricted to potentials of a finite range and can have a wider range 
of applicability. The Krein-Friedel formula may be used even for long-ranged potentials 
when the sum over phase shifts is properly regularized. In the case of the Aharonov-Bohm 
potential it is the (^-function regularization that gives the correct answer. By means of the 
Krein-Friedel formula the change of the DOS induced by the Aharonov-Bohm potential has 
been calculated for different self-adjoint extensions which correspond to different physics 
inside the flux tube. For the conventional setup, i. e., with zero boundary conditions at 



the boundary of the flux tube, our result for the DOS p^ ) conflrms the anticipation of 
Comtet, Georgelin, and Ouvry that the change of the DOS is concentrated at the zero 
energy. Whenever a bound state is present it is always accompanied by a resonance at 
an energy proportional to the absolute value of the energy of the bound state. In the 
presence of a bound state or a phase-shift flip, the symmetry of the spectrum with respect 
to the change a — a ± 1 is broken. Nevertheless, scattering cross sections continue to be 
periodic with respect to the substitution a — > a ± 1. 

In the case of a regular flux tube of a flnite radius i?, the question of the number of 
bound states has been clarifled. It has been shown that the number depends not only 
on the total flux but also on the energy of magnetic fleld. It means that one has to be 
careful in the choice of the regularization when discussing the physics with anomalous 
magnetic moment (^m > 2 or with a short-range interaction inside the flux tube. When 
the AB vortex arises by some physical mechanism such as vortices in the superconductor 
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of type II do, then the profile of a magnetic field is determined by the condition of the 
minimum of energy. Our considerations have only been quantitative. The open question 
which remains is what is the qualitative dependence on the energy and whether (and if 
then, how) the number of bound states depends on higher moments / B^{r)d?Y of the 
magnetic field B{r), where n > 3. We have shown that a bound state in the / = —n — 1 
channel can exist provided l — rj is sufficiently small. In the discussion of the R ^ limit, 
the existence of a critical coupling = "2 was established. Provided the coupling with the 
interior of the flux tube is smaller or greater than the critical one, or when the coupling 
is not renormalized, the limit coincides with that of the impenetrable flux tube. Thereby, 
the result of Ref. HS] for = is generalized. The Aharonov-Casher theorem about 
the number of zero modes has been corrected. It has been shown that they merge with 
the continuous spectrum in the i? — > limit. The origin of the phase-shift flip has been 
investigated. It was shown that the phase-shift flip [O] may occur even in the absence of 



spin. The conditions for the phase-shift flip formulated in Refs. [14, 53] have been found 



necessary but not sufficient. In particular, provided no bound state is in the / = —n — 1 
channel for R ^ 0, then the phase-shift flip will not occur at the R ^ limit. 

A quantum-mechanical and nonrelativistic criterion of stability ( |114| ) in 2 + 1 dimen- 
sions has been discussed that only involves fundamental parameters of matter. Despite 
that this criterion was derived under very restrictive hypotheses, it was shown to be com- 
patible with other, quantum- field theory treatments An open question remains 
regarding what will happen if the dynamics of the shrinking of the radius of a flux tube 
is included. Another open question is what fleld configuration actually optimalizes the 
right-hand side of Eq. (|114|) , i. e., what is the infimum of the right-hand side of Eq. (|114|) 
with respect to variations of magnetic field B{r), subject to the constraint (|89|) . 

It has been shown that in the presence of bound states or a phase-shift fiip, the 
differential scattering cross section becomes asymmetric with respect to the substitution 
(p —>■ —(p, and gives rise to the Hall effect. The Hall resistivity has been calculated [see 
Eq. ( |118| )]. The persistent current of free electrons in the plane has been discussed and 
the results of Ref. ^ have been extended. In contrast to Ref. we have found that the 
total persistent current in the spin one-half case is different from zero, and the persistent 
current of both the spinless and the spin one-half fermions depends linearly on rj. The 



49 



formal similarity between the scattering of electrons in the AB potential and that in the 
field of a cosmic string suggests that a persistent current will also appear in the latter 
case. 

The known results for the 2nd virial coefficients of the nonrelativistic anyons 0, ^ 



were reproduced and their values were calculated for the case when anyons interact via a 
short-range attractive potential, f/(r — r'), proportional to the Dirac 5 function given by 
Eq. ( |131| ) . The case considered here is complementary to the case of a long-range potential 
g/r'^i considered by Loss and Fu |^2[. In the present case, the 2nd virial coefficients were 



shown to be remarkably stable when such an interaction is switched-on. They do not 
change when the interaction is switched on until the coupling constant reaches its critical 
strength. At the critical coupling the 2nd virial coefficients become discontinuous as a 
function of a but their periodicity with respect to a — > a ± 2 and the interpolation 

{a± 1) = a2^{a) between different statistics still hold. These results were obtained by 
using the ^-function regularization which has an obvious advantage with regard to other 
finite box or harmonic potential regularizations. By using the (^-function regularization 
one avoids the necessity of a discretization of energy levels for computing the partition 
function. Our results naturally generalize those of Blum et al. for nonrelativistic spin 
one-half anyons for the case of an anomalous magnetic moment and when a bound state 
occurs in the spectrum. 

Our results have been formulated in terms of self-adjoint extensions. A self-adjoint 
extension is the rigorous limit i? — 0. Therefore, in the case of a flux tube of flnite radius 
i?, our results can be applied when the radius R is negligibly small when compared to all 
other length scales in the system under consideration. The parameters A_„ and A^„_i of 
the relevant self-adjoint extension are then determined from bound state energies in the 
/ = —n and / = —n — 1 channels. Hence, to check some of our results experimentally it is 
not necessary to use a singular magnetic fleld. By using the duality discussed in Sec. |^the 
physics in the presence of an almost singular flux tube can be tested in experiments with 
slow electrons scattered off a general rotationally invariant two-dimensional magnetic fleld 
B{x, y) that obeys the flnite-flux condition (PU]). From the experimental point of view we 
stress the following experiments that should be performed: 

• Measurement of the phase shift flip in the I = —n — 1 channel. Depending on the 



50 



flux it may or may not occur in this channel. 



• Measurements of the Hall resistivity either in the dilute vortex limit or for the single 
vortex. 

• Measurements of the resonance. 

A hypothetical setup to observe the resonance is to study the transmission through a 
single flux tube. A realistic physical realization is that suggested originally by Rammer 
and Shelankov [jll| and later realized experimentally by Bending, Klitzing, and Ploog 
10 | (see Fig. |l]). When homogeneous magnetic field is switched on, the conventional 



superconductor is penetrated by vortices of flux with a = 1/2, i. e., exactly that flux 
at which the resonance is at infinity. Recent measurements on YBCO-delta rings with 
three grain-boundary Josephson junction p3], however, reported the observation of vor- 



tices that carry a flux a = 1/4 which is sma//er than the standard flux quantum hc/2e 
(corresponding to a = 1/2) in the superconductor. Therefore, when the high-T^ YBCO 
film is used as a gate on top of the heterostructure containing 2DEG, the resonance is at 
some finite energy and could in principle be observed. Another possibility is to measure 
energy dependence of phase shifts and (as a precursor of the resonance) the time delay 

dZI). 

An interesting open problem concerns the shape of the resonance (^3]) in the AB 
scattering that is not of the Breit-Wigner form. The latter is a general consequence of 
the analyticity of scattering amplitudes that is usually taken as an equivalent to causality. 
However, the rigorous proof of the equivalence requires precise localization in time and 
in the energy for the incident wave packet. This is, however, impossible due to the 
uncertainity relations impossible We have undertaken the analysis of the problem in 
terms of the Jost functions. However, we have not found the origin of this behavior and 
we postpone the solution to this problem. The discussion might in principle shed some 
light on the analyticity principle usually adopted in the axiomatic quantum-field theory 
as a substitute for causality. 

Another interesting problem is to find the analog of the Levinson theorem (see Ref. 
4T| , p. 356, and Ref. [^]) and the generalization of Bargmann's inequalities for 



singular potentials. In the case of regular potentials of finite range one knows, thanks to 
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this theorem, that if the phase shifts are normahzed so that 6i{E) = for E oo, then 
5i(0) gives the number of bound states in the channel /. In the case of the AB potential 
the phase shifts (|^) essentially do not depend on the angular momentum / and on the 
energy E. In the case of a < 1 they are either — vro; for positive I or na for negative I. It 
is impossible in general to normalize phase shifts in the singular case in the same way as 
for regular potentials, since they do not depend on the energy. 
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A Some useful integrals 



To perform the integral in Eq. (|14|) one uses the result that 

1 r. 



(157) 



and deforms the integration contour as in Fig. |. and if^^) 

{z) are the Hankel 

functions, 
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sm UTC 
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For real arguments. 
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(158) 



(159) 



[see Ref. [^, relations (9.1.3) and (9.1.4)]. The integral with H^\kry) is now deformed to 
the contour Ci while that with H'^\kry) is deformed to the contour C2. Note that H^\z) 
[H^\z)\ is exponentially decreasing in the upper (lower) half of the complex plane (see 
Appendix |B|). Since it is assumed that < Vy the integrand in Eq. ( p!^ is exponentially 
decreasing as well. Now, the contributions of any of the two integrals from the imaginary 



52 



axis do not vanish. However, by using the identities (9.1.35) and (9.1.39) from Ref. p^ , 

Mze^^')=e^''^'Mz), (160) 

i7(2)(^e-'^^) = -e""i7«(^), (161) 

i7(i)(^e"0 = -e-''''Hl^\z), (162) 

one can show that 

'"^^ Uzr,)H^^\zry) = - f -^^J^(zr,)Hi'\zry), (163) 



loo (/ — Z J ~ioo (/ — Z 

i. e., the integrals on the imaginary axis cancel in the sum. Thanks to ( |161| ) one has for 
z/ = 

H^^\ze-n = -Hl^\z), (164) 



and the same applies to the integral in Eq. ([121) , too. Once the integrals in Eqs. ([T^ ) and 
(|T^) are represented as the sum of two contour integrals in the complex plane they can 
be simply taken by the residue theorem. 

In order to calculate Tr AGq,(x, x, M) note that 



Tr AG'a(x, X, M) = / /\GJ^, x, M) d 



^x 



' ' du^7—,—-—ri^,^—j-.. (165) 



'— oo J —oo (coshcj + cosh'i?)^ 1 + 
To perform the integral here one makes the substitution 



^"^~2^' ^^^^^ 



with a Jacobian equal to 2. After some manipulations one finds 



oo rao 

dd / duo 



oo j-oo (coshcij + cosht?)^ 1 + e^' 
°° dx f°° cosh(2r7 — l)y 



lo cosh^ X Jo cosh^ y 
According to formula (3.512.1) of Ref. |3^ 



dy. (167) 



r^^5!^^^^rf!, = 2B(l + ,.2-„), (168) 

Jo cosh y 
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where B{x,y) is the Euler beta function, 

B{x,y) 



and r(x) is the Euler gamma function Since 1) = 1, 

roo 



cosh X 



dx = 1. 



(169) 



(170) 



Now, r(3) = 2 and 



25(1 +i],2-r]) = T{1 + r])T{2 - r]) = 

rJ{l-l^)T{l^)T{l-r]). (171) 

Therefore, by using ( p.78| ), 

Jo cosh y sin ttt] 

B Cylindrical functions and the exterior solution 

The asymptotic behavior of Bessel functions and their derivatives at infinity as z ^ oo 



can be found directly in Ref. |3^. The asymptotic behavior of the Bessel function Jiy{z) 
and their derivatives at the origin as 2; —> is determined by relation (9.1.10) of Ref. ||34|| . 
The asymptotic behavior of the Hankel function at the origin for |z/| < 1 can be calculated 
from the asymptotic behavior of Ju{z) by using relation (9.1.3) of Ref. 
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(173) 
(174) 



where 7 is Euler's constant. The asymptotic behavior of Kyi^z) for z/ < 1 when 2; ^ is 
determined with the help of the last two formulas and relations (9.6.4) and (9.6.6) of Ref. 



3^ , that hold for — vr < arg z < 7t/2. Therefore, 
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r(i-i/) fzy 



T{l + iy) \2 



+ Oiz''") 



V] z 



r(i + z/) V2 



\0{z 



2-u\ 



(176) 
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+ 0{z 



l-v\ 



otherwise. Here we have used the identity = r(l + u) together with 

TT 



The wave function outside the flux tube is 



Sm UTT 



(177) 



(178) 



(179) 



For its logarithmic derivative rKl[r) / Ky[r) one has in general (see formula 9.6.26 of Ref. 
0) 

(180) 



K{z) 



-V — z -— = 1^ — 2;- 



K,{z) K,{z) K,{z) 

Its asymptotic behavior as z ^ depends whether u is less than, equal to, or greater 

than 1. One has 

K'Jz) „ r(l-z/) fz^^^ 



for < z/ < 1/2, 
for l/2<u < 1, 
for z/ = 1, and 



K'iz) 



-v-2v 



T{l + u) V2 



+ 0{z 



KJz) 



'KJz) 



for V > 1. At infinity as r ^ oo 



Til-v) fz\'^'' 
~-l+p^z'lnz + 0(.-') 



2(^-1) 



z' + 0{z'\ z 



KJr) 



(181) 
(182) 

(183) 

(184) 

(185) 



C Hypergeometric function and the interior solution 

In the case of homogeneous field regularization, the wave function inside the fiux tube of 



radius R at the energy E is given by (cf. Ref. p. 458) 

V'(r) 



M 



e-€/2^|i|/2 
f\l\+l e(r) 



1 



2;2 



(186) 
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Here ^ is the flux within the radius r in units of the flux quantum $0; ^ = ^s'^^/S = 
<l>(r)/$0; Ib = {hc/\e\By^'^ is the magnetic length, fc^ = 2mE/h^, and the magnetic 
moment coupling has been assumed inside the flux tube. Since we are interested in bound 



states, the parameter x = kR will be introduced in analogy with Ref. where, as above, 
K = \j2m\Ei,\/?i^ , Eh being the bound state energy. One can show that 

„2 



X 

2a 



2mEh 



ft 



— '■B 



In. 
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The logarithmic derivative of (|186|) at the flux tube radius r = i? is then 



-a + \l\ + 2a 



r=R 



a M {a + l,b + l,a) 
b M{a,b,a) ' 



where a and b are given by Eq. { ^7\ 



The Kummer function M(a, 6, z), frequently denoted by iFi(a, b, z) [34|, is defined by 



M(a.M) = l + -. + A__2_ + ... + + 



where 



;i89) 



a„ = a(a + 1) . . . (a + n — 1), 
bn = b{b+l)...{b + n-l). 



(190) 



Due to the formula 13.4.8 of Ref. 



-^M{a, b, z) = ^M{a + l,b + l,z) 
dz b 



Using the simple fact that, provided 6 > 0, 

a < b =^ 

one finds for real z = c > 



a a + 1 

b ^ b+r 



(191) 



(192) 



M(a, 6, c) < M(a + 1, 6 + 1, c). 



(193) 



Note that in the cases discussed here [see Eqs. (0) and (|87D] one has not only 6 > 
but even 6 > 1. Moreover, provided there is neither magnetic moment coupling nor an 
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attractive potential inside the flux tube one also has a > [see Eq. (|81|)], and Eq. ( ^2]) 
holds. However, one cannot maintain 

M{a,b,c)>0 (194) 

when an attractive potential of whatever origin is inside the flux tube. For example, for 
gm — "2 = 2e < 4 and x G 0, ^/2ae^ [see Eq. (W^] one has — 1 < a < 0. In this case 
a + 1 > and M{a + 1, 6 + 1, c) > 0. From the definition (|189|) of M(a, 6, c) one finds 
easily for such a that 

M{-\a\,b,c) = 2- M{\a\,b,c). (195) 

Therefore, if the relation 

M{\a\,b,c)>2 (196) 



is satisfied, one finds M(a, b, c) < 0. Hence, 

M(a + 1,6+ l,c) 



< 0. (197) 



M(a,6, c) 

The ratio at a; = is nothing but a parameter ai of Ref. [^. It is here where the 
error |^ in Ref. |^ is made since they claimed ai to be positive for whatever are the 
parameters a, b, and c. For our purposes it is more important to discuss the property of 
the ratio given in Eq. (^2)) in place of ( |197| ). Note that the ratio stays positive, 

2aM(a + l.t+l.a) 

b M{a,b,a) ^ ' 

provided Eq. ( |196| ) holds. Therefore, at those values of x [it enters Eq. ( |196| ) via a] it is 
impossible to satisfy Eq. (Q) with the homogeneous field regularization. 

D C, functions 

The Riemann C, function is defined by 

oo 

Us)--=Y.^~'- (199) 
1=1 

It is an analytic function with a simple pole at s = 1. The Hurwitz C,h{SiX) function, 
a; ^ Z, is a generalization of Cr{s), 

oo 

Ch{s,x) :=^(/ + x)-^ (200) 

1=0 
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For s = — 1 



and 

C,^(-l, x) = lx{l -x)-^ = - ^) + U-l)- (202) 
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